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I Introduction 

The theory of multiphase polymer systems has a venerable tradition. The 'classical' 
theory of polymer demixing, the Flory-Huggins theory, was developed already in the 
forties of the last century.'^'l^ It is still the starting point for most current approaches 
- be they improved theories for polymer (im)miscibility that take into account the 
microscopic structure of blends more accurately, or sophisticated field theories that 
allow to study inhomogeneous multicomponent systems of polymers with arbitrary 
architectures in arbitrary geometries. In contrast, simulations of multiphase polymer 
systems are relatively young. They are still limited by the fact that one must simulate 
a large number of large molecules in order to obtain meaningful results. Both powerful 
computers and smart modeling and simulation approaches are necessary to overcome 
this problem. 

In the limited space of this chapter, I can only give a taste of the state-of-the-art in 
both areas, theory and simulation. Since the theory has reached a fairly mature stage 
by now, many aspects of it are covered in textbooks on polymer physics. '^^'^'^ The 
information on the state-of-the art of simulations is much more scattered. This is why 
I have put some effort into putting together a representative list of references in this 
area - which is of course still far from complete. 

The chapter is organized as follows. In Section |TI1 I briefly introduce some basic 
concepts of polymer theory. The purpose of this part is to make the chapter accessible 
to readers who are not very familiar with polymer physics; it can safely be skipped 
by the others. Section IIIII is devoted to the theory of multiphase polymer systems. I 
recapitulate the Flory-Huggins theory and introduce in particular the concept of the 
Flory interaction parameter (the x parameter), which is a central ingredient in most 
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theoretical descriptions of multicomponent polymer systems. Then I focus on one of 
the most successful mean-field theories for inhomogeneous (co)polymer blends, the self- 
consistent field theory. I sketch the main idea, discuss various aspects of the theory 
and finally derive popular analytical approximations for weakly and strongly segregated 
blends (the random phase approximation and the strong-segregation theory). 

In Section IIVI I turn to discussing simulations of multiphase polymer systems. 
A central concept in this research area is 'multiscale modeling': Polymers cannot be 
treated at all levels of detail simultaneously, hence coarse-grained models are used in 
order to study different aspects of the materials in different simulations. This allows 
one to push the simulation limits to larger length and time scales. I describe some of 
the most popular coarse-grained structural and dynamical models and give an overview 
over the state-of-the-art of simulations of polymer blends and copolymer melts. 

II Basic Concepts of Polymer Theory 

For the sake of readers who are not familiar with polymer theory, I begin with recapit- 
ulating very briefly some basic concepts. 

Polymers are macromolecules containing up to hundreds of thousands of atoms. At 
first sight, one would not expect such molecules to be easily amenable to theoretical 
modeling; however, it turns out that the large size of the molecules and their highly 
repetitive structure in fact simplifies things considerably. Since polymer molecules 
interact with many others, details of local interactions average out and polymers can 
often be characterized by a few effective quantities, such as their topology, the local 
stiffness along the backbone, the bulkiness, the compatibility/incompatibility of the 
building blocks etc. Already decades ago, pioneers like Flory,'^ Edwards,'^ de Genne^ 
have established theoretical polymer science as a highly successful field of research, 
which brings together scientists from theoretical chemistry, statistical physics, materials 
science, and even the biosciences, has created a wealth of new beautiful theoretical 
concepts, and has not lost any of its fascination for theorists up to date. 

II. 1 Fundamental Properties of Polymer Molecules 

The characterizing property of polymers is their highly modular structure. They are 
composed of a large number of small building blocks (monomers), which are often all 
alike, but may also be combined to arbitrary sequences (in the case copolymers and 
biopolymers). The monomers are arranged in chains, which are usually flexible on 
the nanometer length scale, i.e., they can form kinks at little energetic expense, they 
curve around and may assume a large number of conformations at room temperature. 
The properties of such flexible polymers are largely determined by the entropy of the 
chain conformations. For example, the number of available conformations is reduced if 
molecules are stretched, which leads to a purely entropic restoring spring force^*^ (rubber 
elasticity). Exposed to stress, polymeric systems respond by molecular rearrangements, 
which takes time and results in time-dependent strain (viscoelasticity) . 
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The fundamental processes that govern the behavior of polymeric materials do not 
depend on the chemical details of the monomer structure. For qualitative purposes, 
polymer molecules can be characterized by a few properties such as 

• The architecture of the molecules (linear chains, rings, stars, etc.) 

• Physical properties (local chain stiffness, chain size, monomer volume) 

• Physicochemical properties (monomer sequence, compatibility, charges) 

• Special properties {e.g., a propensity to develop crystalline or liquid crystalline 
order) . 

II. 2 Coarse-Graining, Part I 

The notion of 'coarse-graining' has lately become a buzzword in materials science, 
but the underlying concept is actually quite old in polymer science. The need for 
coarse-graining results from the fact that polymeric materials exhibit structure on very 
different length scales, ranging from Angstrom (the monomeric scale) to hundreds of 
nanometers (typical molecule extensions) or micrometers (supramolecular aggregates). 
It is not possible to treat all of them within one common theoretical framework. There- 
fore, different theoretical descriptions have been developed that deal with phenomena 
on different length and time scales. On the microscale, chemical details are taken 
into account and the polymers are treated at an atomistic level. This is the realm of 
theoretical chemistry. On the mesoscale, simplified molecule models come into play 
(string models, lattice models, bead-spring models, see below), whose behavior can be 
understood with concepts from statistical physics. Finally, on the macroscale, poly- 
meric materials are described by continuous fields (composition, strain, stress etc.) 
with certain mechanical properties, and their behavior can be calculated with methods 
borrowed from the engineers. 

In the following, I shall mainly focus on the mesoscale level, where polymers are de- 
scribed by extended molecules made of simplified "monomeric" units, each representing 
several real monomers. Even within that level, one still has some freedom regarding 
the choice of the coarse-grained units. This is illustrated in Fig. 13. H where different 
coarse-grained representations of a polymer are superimposed onto each other. Poly- 
mers have remarkable universality properties, which allow one to link different coarse- 
grained representations in a rather well-defined manner, as long as the length scales 
under consideration are much larger than the (chemical) monomer length scale. For 
example, starting from one (atomistic or coarse-grained) model, we can construct a 
coarse-grained model by combining m "old" units to one "new" unit. If m is suffi- 
ciently large, the average squared distance (cP) between two adjacent new units will 
depend on m according to a characteristic power law 
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Figure 3.1: Mesoscopic coarse-grained representations of a polymer molecule. Light 
pearl necklace in the background: Bead chain with bonds of fixed length. Solid and 
dashed lines, chain of coarse-grained units linked by bonds of variable length. 

where the exponent v depends on the environment of a polymer, but not on chemical 
details.'^'^ In a dense polymer melt, one has u = 1/2 (see below). Similar scaling laws 
can be established for other chain parameters. 

II.3 Ideal Chains 

In mesoscopic polymer theories, one often uses as a starting point a virtual polymer 
chain where monomers that are well separated along the polymer backbone do not 
interact with each other, even if their spatial distance is small. Such polymers are 
called 'ideal chains'. Even though they are mere theoretical constructions, they provide 
a good approximative description of polymers in melts and in certain solvents ('Theta'- 
solvents, see below). 

II. 3.1 A Paradigm of Polymer Theory: The Gaussian Chain 

Let us now consider a flexible ideal chain with monomers, which we coarse- grain 
several times as sketched in Fig. 13. 1|, until one coarse-grained monomer unites m 'real' 
monomers. For large m, the resulting chain is a random walk in space consisting 
of uncorrelated random steps dj of varying length. According to the central limit 
theorem of probability theory,^' the steps are approximately Gaussian distributed, 
'P(d) ~ exp(— (i^/2m(T^), where a does not depend on m. The same random walk 
statistics can be reproduced by a Boltzmann distribution with an effective coarse- 
grained Hamiltonian 

N/m 

'^m = \^,Y.^l (3-2) 

i=l 
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The Hamiltonian Ti^ describes the energy of a chain of springs with spring constant 
kgTm/a^. The coarse-graining procedure has thus ehminated the information on chemi- 
cal details (they are now incorporated in the single parameter a), and instead unearthed 
the entropically induced elastic behavior of the chain which lies at the heart of rubber 
elasticity. Eq. ()3.2p is also an example for universal behavior in a polymer system (see 
section |IL21): The coarse-grained chain is self-similar. Every choice of m produces an 
equivalent model, provided the spring constant is rescaled accordingly. The distance 
between two coarse-grained units exhibits a scaling law of the form ()3.ip as a function 
of m, {(f) = Sa^m^" with u = 1/2. 

Based on these considerations, it seems natural to define a 'generic' ideal chain 
model based on Eq. ()3.2p with m = 1, 



N-1 



Y,i^^+l-^^)^ (3-3) 



the so-called 'discrete Gaussian chain' model. For theoretical purposes, it is often 
convenient to take the continuum limit: The index i in Eq. (j3.2p . which counts the 
monomers along the chain backbone, is replaced by a continuous variable s, the chain 
is parametrized by a continuous path R(s), and the steps d correspond to the local 
derivatives dR/ds of this path. The effective Hamiltonian then reads 



This defines the continuous Gaussian chain. The only material parameters in Eq. ()3.4p 
are the chain length and the 'statistical segment length' or 'Kuhn length' b. Even 
those two are not independent, since they both depend on the definition of the monomer 
unit. An equivalent chain model can be obtained by rescaling N — )■ N/X and 6^ — )• b^X. 
Hence the only true independent parameter is the extension of the chain, which can be 
characterized by the squared gyration radius 



1 



N 



Rl = j^i ds (R(s) - R)2 = h'N/Q, (3.5) 







where R = J ds R(s) is the center of mass of the chain. The quantity Rg sets 
the (only) characteristic length scale of the Gaussian chain, and all length-dependent 
quantities scale with Rg. For example, the structure factor is given by 



^(k) = ^( 



N 







Ng^ik'Rl), (3.6) 



with the Debye function 

9^{x) = ^{e'^ -l + x). (3.7) 

The Gaussian chain is not only a prototype model for ideal chains, it also provides 
a general framework for mesoscopic theories of polymer systems. The Hamiltonian 
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Tig (Eq. (j3.4p ) is then supplemented by additional terms that account for interactions, 
external fields, constraints {e.g., chemical crosslinks) etc. In this more general context, 
the Hamiltonian (|3.4p is often referred to as 'Edwards Hamiltonian'. 

Finally in this section, let us note that from a mathematical point of view, the 
probability distribution of chain conformations defined by Eq. (j3.4p is a Wiener mea- 
sure.'^'^ The continuum limit leading to Eq. (13. 4p is far from trivial, but well-defined. 
I shall not dwell further into this matter. 

II.3.2 Other Chain Models 

The Gaussian chain model is a common starting point for analytical theories of long 
flexible polymers on sufficiently large length scales. On smaller length scales, or for 
stiffer polymers, or for computer simulation purposes, other types of coarse-grained 
models have proven useful. I briefly summarize some popular examples. 

The wormlike chain model is a continuous model designed to describe stiff polymers. 
They are represented by smooth paths R(s) with flxed contour length N, where 
the parameter s runs over the arc length of the curve, i.e., the derivative vector 
u = dR/ds has length unity, |u| = 1. The paths have a bending stiffness rj, such 
that they are distributed according to the effective Hamiltonian 



The wormlike chain model is particularly useful if local orientational degrees of 
freedom are important. 

The freely jointed chain is a discrete chain model where the chain is composed of N 
links of flxed length. It is often used to study general properties of ideal chains. 

The spring-bead chain is a chain of beads connected with springs. It has some resem- 
blance with the discrete Gaussian chain model, except that the springs have a 
flnite equilibrium length. Spring-bead models are popular in computer simula- 
tions. 

In lattice models, the monomer positions are confined to the sites of a lattice. This 
simplifies both theoretical considerations and computer simulations. 

II. 4 Interacting Chains 

The statistical properties of chains change fundamentally if monomers interact with 
each other. Such interactions are readily introduced in the coarse-grained models pre- 
sented above. In the discrete models, one simply adds explicit interactions between 
monomers. In the continuous path models, one supplements the energy contribution 
for individual ideal chains, Eq. (13. 4p or (13. Sp . by an interaction term, such as 






(3.9) 
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(for weak interactions), where the 'monomer density' p{r) is defined as 

p(r) = V/ ds6ir-R^is)) (3.10) 
a -^0 

and the sum a runs over the polymers Rq:(s) in the system. Eq. (j3.9p corresponds to a 
virial expansion of the local interaction energy in powers of the density. In many cases, 
only the quadratic term (v) needs to be taken into account ('two-parameter Edwards 
model'). The Ansatz (j3.9p is suitable for dilute polymer systems - dense systems 
are discussed below (Sec. III.4.2p . The generalization to multiphase systems where 
monomers may have different type A,B,- • • is straightforward. One simply operates 
with different densities PAiPb-,' ' ' a-nd interaction parameters vaa-, vab, ■ ■ ■ • 

Interactions complicate the theoretical treatment considerably and in general, exact 
analytical solutions are no longer available. The properties of interacting polymer sys- 
tems have been explored theoretically within mean-field approximations, renormalization- 
group calculations, scaling arguments, and computer simulations. To set the stage for 
the discussion of multiphase systems in sections IIIII and IIVI I will now briefly sketch 
the most important scenarios for monophase polymer systems. 

II. 4.1 Polymers in Solution and Blobs 

We first consider single, isolated polymer chains in solution. Their properties depend 
on the quality of the solvent, which is incorporated in the second virial parameter v in 
Eq. (|3.9p (the three-body parameter w is typically positiv^^^. In good solvent (f > 0), 
monomers effectively repel each other, and the chain swells. Extensive theoretical 
worlPlhas shown that the scaling behavior (Eq. ()3.ip ) remains valid, but the exponent 
u increases from = 1/2 for ideal chains to ~ 3/(d + 2), where d is the spatial 
dimension (more precisely, = 0.588 in three dimensions). This is the famous 'Flory 
exponent', which characterizes the scaling behavior of so-called 'self-avoiding chains'. 
Accordingly, the gyration radius of the chain scales with the chain length like 

Rl ~ N^". (3.11) 

In poor solvent {v < 0), monomers effectively attract each other and the chain collapses. 
At the transition between the two regimes, the 'Theta point' (vq ~ 0), the scaling 
behavior basically corresponds to that of ideal chains (z^ = 1/2), except for subtle 
corrections due to the three-body w-teim^ 

Eq. ()3.1ip describes the behavior of single, unperturbed chains. Even in good sol- 
vent, the self-avoiding scaling is often disturbed. For example, the chains cannot swell 
freely if they are confined, or if they are subject to external forces. Another important 
factor is the concentration of chains in the solution: If many chains overlap, the in- 
trachain interactions are screened on large length scales. Loosely speaking, monomers 
cannot distinguish between interactions with monomers from the same chain and from 
other chains. As a result, chains no longer swell and ideal chain behavior is recovered. 

This mechanism applies in three or more spatial dimensions. Two dimensional chains 
segregate.ElEIISl 
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a) b) C) 




Figure 3.2: Illustration of the blob model in different situations: a) concentrated poly- 
mer solution b) Polymer confined in a slit c) chain confined in a spherical cavity d) 
structure formation in solutions of miktoarm star copolymers (after Ref . |16p . See text 
for explanation. 

All of these situations can be analyzed within one single ingenious framework, the 
'blob' picture introduced Daoud et al. in 1975.'^ It is based on the assumption that 
there exists a crossover length scale ^ below which the chain is unperturbed. Blobs are 
volume elements of size ^ within which the polymers behave like self-avoiding chains. 
On larger scales, the polymer behaves like an ideal chain consisting of a string of blobs. 
Every blob contains m ~ ^^^'^ monomers and carries a free energy of the order k^T. 
These simple rules are the whole essence of the blob model. I shall illustrate their use 
by applying them to a number of prototype situations depicted in Fig. 13. 2i 

Concentrated polymer solution (Fig. l3.2l a). For polymer concentrations (p, we calculate 
the crossover length scale ^ from self-avoiding to ideal behavior. Since ^ is the 
blob size, we can simply equate 4> = m/^^, i.e., ^ ~ 

Polymer confined in a slit (Fig. 13.21 b) . We consider the free energy penalty F on the 
confinement. Here, the blob size is set by the width R of the slit. Each blob 
contains m ~ R^^'^ monomers, hence the total free energy scales like F ^ N/m ~ 
NR-^I''. 

Polymer confined in a cavity (Fig. 13.21 c) The result b) also holds for chains con- 
fined in a tube. In closed cavities, however, the situation is different due to 
the fact that the cavity constrains the monomer concentration. The resulting 
blob size is i ~ {N/R^Y/^^-'^''\ and the free energy of confinement scales as 
F ~ (i?/^)^ ~ (R/N'^)^/^^~'^'^\ This has been discussed controversially, but was 
recently confirmed by careful computer simulations.'^ 

ABC miktoarm star copolymers in selective solvent (Fig. 13.21 d). Last, I cite a recent 
application to a multiphase polymer system. Zhulina and Borisov^*' have stud- 
ied ABC star copolymers by means of scaling arguments. They derived a rich 
state diagram, according to which ABC star copolymers may assemble to several 
types of nanostructures, among other spherical micelles, dumbbell micelles, and 
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striped cylindrical micelles. This is only one of numerous examples where scaling 
arguments have been used to analyze complex multicomponent systems. 

II.4.2 Dense Melts 

Dense melts can be considered as extreme cases of a very concentrated polymer solu- 
tion, hence it is not surprising that the chains effectively exhibit ideal chain behavior. 
In fact, the situation is more complicated than this simple argument suggests. The 
quasi-ideal behavior results from a cancellation of two effects: On the one hand, the 
intrachain interactions promote chain swelling, but on the other hand, the chain pushes 
other chains aside ('correlation hole'), which in turn exert pressure and squeeze it. De- 
viations from true ideal behavior can be observed, e.g., at the level of chain orientational 
correlations.'^'^ Nevertheless, the ideality assumption is a good working hypothesis 
in dense melts and shall also be used here in the following. 

II. 5 Chain Dynamics 

In this article I focus on equilibrium and static properties of polymer systems. I can 
only touch on the possible dynamical behavior, which is even more diverse. 

In the time scales of interest, the motion of polymers is diffusive, i.e., the inertia 
of the macromolecules is not important. Three prominent types of dynamical behavior 
have been established. 

In the Rouse regime, the chain dynamics is mainly driven by direct intrachain inter- 
actions. This regime is encountered for short chains. The dynamical properties 
of ideal chains can be calculated exactly, and the results can be generalized to 
self-avoiding chains using scaling arguments. One of the important properties of 
Rouse chains is that their sedimentation mobility does not depend on the chain 
length A^. Hence the diffusion constant scales like D ~ ^/N, and the longest 
internal relaxation time, which can be estimated as the time in which the chain 
diffuses a distance Rg, scales like r ~ N'^'^^^. 

In the Zimm regime, the dynamics is governed by long-range hydrodynamic inter- 
actions between monomers. This regime develops for sufficiently long chains 
in dilute solution. They diffuse like Stokes spheres with the diffusion constant 
D ~ ^/Rg, and the longest relaxation time scales like t ^ Rg. In concentrated 
solutions, the hydrodynamic interactions are screenecP and Rouse behavior is 
recovered after an initial Zimm period.'^ 

The reptation regime is encountered in dense systems of chains with very high molec- 
ular weight. In this case, the chain motion is topologically constrained by the 
surrounding polymer network, and they are effectively confined to move along a 
tube in a a snake-like fashion.SlISI xhe diffusion constant of linear chains scales 
like D ~ and the longest relaxation time like r ~ N^. 
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Spheres (S) Cylinders (C) Gyroid (G) Lamellae (L) 

Figure 3.3: Self-assembled copolymer mesophases. 



This description is very schematic and oversimplifies the situation even for fluids of 
linear polymers. Moreover, most polymer materials are not in a pure fluid state. They 
are often cooled down below the glass transition, or they partly crystallize - in both 
cases, the dynamics is frozen. Chemical or physical crosslinks constrain the motion of 
the chains and impart solid-like behavior. In multiphase polymer systems, the situation 
is further complicated by the fact that the glass point or the crystallization temperature 
of the different components may differ, such that one component freezes where the other 
still remains fluid. '23^22 xhe following discussion shall be limited to fluid multiphase 
polymer systems. 



Ill Theory of Multiphase Polymer Mixtures 

After this general overview, I turn to the discussion of polymer blends. We consider 
dense mixtures, where the polymers are in the melt regime (Sec. III.4.2j) . Moreover, we 
assume incompressibility - the characteristic length scales of density fluctuations are 
taken to be much smaller than the length scales of interest here. 

Monomers of different type are usually slightly incompatible (see Section fllLl-Sp . In 
polymers, the incompatibilities are amplified, such that macromolecules of different type 
tend to be immiscible: Blended together, they demix and develop an inhomogeneous 
multiphase structure where microdroplets of one phase are finely dispersed in another 
phase. 

In order to overcome or at least control this situation, copolymer molecules can be 
added in which the two incompatible components are chemically linked to each other. 
They act as compatibilizer, i.e., they shift the demixing transition and reduce the inter- 
facial tension between different phases in the demixed region. At high concentrations, 
they are found to self-organize into a variety of ordered mesophases (microphase sepa- 
ration; see, e.g., structures shown in Fig. 13. Sp . Hence copolymers can also be used to 
manufacture nanostructured materials in a controlled way. 

Nowadays, the theory of structure formation in polymer blends has reached a highly 
advanced level and theoretical calculations have predictive power, e.g., with respect to 
structures that can be expected in new polymeric materials. In this section, I shall 
present some of the most successful theoretical approaches. 
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III.l Flory Huggins Theory 

I begin with sketching the Flory-Huggins theory, which is the classical theory of phase 
separation in polymer blends, and which in some sense lays the foundations for all later, 
more sophisticated theories of polymer mixtures. 

III. 1.1 Basic Model for Binary Blends 

We consider a binary blend of homopolymers A and B with length Nj^ and Nb, and 
volume fractions and ^b- According to FlorjC and Huggins,'^ the free energy per 
monomer is approximately given by 

^ = §i ln{^A) + p- H<^b) + X^A^B (3.12) 
kgT Na Nb 

with ^A + = 1- The first two terms account for the mixing entropy of the two 
components, and the last term for the (in) compatibility of the monomers. The param- 
eter X is the famous 'Flory Huggins parameter', which will be discussed in more detail 
below. The generalization of this expression to ternary ABC homopolymer blends etc. 
is straighforward, one only needs to introduce several x-pa^'a^eters XAB, Xbc-, and 
XAC- Here we will only discuss binary systems. 

By minimizing the free energy, Eq. ()3.12p . one easily identifies the region in phase 
space where the mixture phase separates into an A-rich phase and a B-rich phase. At 
low values of X) the blend remains homogeneous. Demixing sets in at a critical value 
2Xc = (l/\/iVj+ l/\/iVB)^ for the critical composition ^a,c = 1/(1 + ^/Na/Nb). The 
region of stability of the homogeneous (mixed) blend is delimited by the " binodal" line 
(see Fig. l3.4p . Beyond the binodal, the homogeneous blend may still remain metastable. 
It becomes unstable at the "spinodal", which is defined as the line where the second 
derivative of fp^ in Eq. (j3.12p with respect to ^a vanishes. An example of a phase 
diagram with a binodal and a spinodal is shown in Fig. 13.41 (left). Fig. 13.41 (right) 
demonstrates the shift of the binodal with varying chain length ratio Na/Nb- 

The Flory-Huggins free energy, Eq. (|3.12p . was originally derived based on a lattice 
model, but it can also be applied to off-lattice systems. It does, however, rely on three 
critical assumptions: 

• The polymer conformations are taken to be those of ideal chains, independent of 
the composition (ideality assumption, cf. Sec. 111.4.^ . 

• The melt is taken to be incompressible, and monomers A and B occupy equal 
volumes. 

• Local composition fluctuations are neglected (mean-fleld assumption). 

In reality, none of these assumptions is strictly valid. The polymer conformations do 
depend on the composition, most notably for chains of the minority component. The 
incompressibility assumption is reasonable, but the volumes per monomer are not equal. 
As a consequence, the x-parameter is not a flxed parameter (at fixed temperature), but 
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Figure 3.4: Left: Phase diagram for a binary AB polymer blend with B-chains twice 
as long as A-chains according to the Flory-Huggins theory. Thick solid line shows 
the binodal line (i.e., the demixing line), thin dashed line the spinodal line {i.e., the 
line where the homogeneous blend becomes unstable). Right: Binodals for binary AB 
blends with different chain length ratios as indicated. 



depends on the composition of the blend (see Sec. IIII.1.3p . Finally, the composition 
fluctuations shift phase boundaries and may even fundamentally change the phase 
behavior, (see Sec. IIII.231) . 



III. 1.2 Inhomogeneous Systems: Flory-Huggins-de Gennes Theory 

Eq. (I3.12p . only describes homogeneous systems. The simplest approach to generalizing 
the Flory-Huggins theory to inhomogeneous systems, e.g., polymer blends containing 
interfaces, consists in adding a penalty on composition variations (V^a)'^ = (V<^b)^- 
The coefficient of the square gradient term can be derived within a more advanced 
mean-field treatment, the random phase approximation, which will be described further 
below (Sec. IIII.3.ip . One obtains the Flory-Huggins-de Gennes free energy functional 
for polymer blends, 

T,^.o[^A{r)] =Poj dr |/,,(cl>^(r)) + + ^)(V<J>a)^} (3.13) 

(with $5 = 1 — ^a)i where 6^ and bs are the Kuhn lengths of the homopolymers A 
and B, and /^^(<I>yi)) is given by Eq. (j3.12p . The functional (j3.13p can be applied if 
composition variations are weak, and have characteristic length scales of the order of 
the gyration radius of the chains ('weak segregation regime', see Sec. IIII.3T]) . 

A very similar functional can be derived in the opposite case, where A- and B- poly- 
mers are fully demixed and separated by narrow interfaces. In this 'strong segregation' 
regime, the blend can be described by the functional (see Sec. IIII.3.2P 

J-3sJ$A(r)] = po / dr|x<I>A$B + ^(^ + ^)(V$^)2|. (3.14) 

At strong segregation, the mixing entropy terms in /fh, Eq. (I3.12p . can be neglected, 
hence the functionals (j3.13p and (j3.14p are identical except for the numerical prefactor of 
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the square gradient term. In the strong segregation hmit, the square gradient penalty 
results from an entropic penalty on A and B segments due to the presence of the 
interface, whereas in the weak segregation limit, it is caused by the deformation of 
whole chains. 

III. 1.3 Connection to Reality: The Flory-Huggins Parameter 

In the Flory-Huggins theory, the microscopic features of the blend are incorporated in 
the single Flory-Huggins parameter x- Not surprisingly, this parameter is very hard to 
access from first principles. 

In the original Flory-Huggins lattice model, x is derived from the energetic inter- 
actions between monomers that are neighbors on the lattice. The interaction energy 
between monomers i and j is taken to be characterized by energy parameters e^j. The 
X- parameter is then given by 



where z is the coordination number of the lattice. It is reduced by two {z — 2) in order 
to account for the fact that interactions between neighbor monomers on the chain are 
fixed and have no influence on the demixing behavior. 

In reality, the situation is not as simple. The miscibility patterns in real blends 
tend to deviate dramatically from that predicted by the Flory-Huggins model. Several 
blends exhibit a lower critical point instead of an upper critical point, indicating that 
the demixing is driven by entropy rather than enthalpy. The critical temperature Tc 
of demixing often does not scale linearly with the chain length N, as one would expect 
from Eq. (13.15P (see Fig. 13. 4p . In blends of polystyrene and poly (vinyl methyl ether) 
(PS/PVME), for example, Tc is nearly independent of N, and the critical concentration 
$c is highly asymmetric even for roughly equal chain lengths, in apparent contrast with 



Formally, this problem can be resolved by arguing that the Flory-Huggins expression 
for X, Eq. (j3.15p . is oversimplified and that x^b^ is not a constant. Several factors 
contribute to the ^-parameter, leading to a complex dependence on the temperature, 
the blend composition, and even the chain length. 

Monomer incompatibility : Monomers may be incompatible both for enthalpic and 
entropic reasons. For example, consider two nonpolar monomers i and j. The 
van-der Waals attraction between them is proportional to the product of their 
polarizabilities a, hence eij oc aiUj and x {o^A — oibY'- More generally, the en- 
thalpic incompatibility of monomers i and j can be estimated by xh oc {5a — ^b)"^, 
where the 5i are the Hildebrand solubility parameters of the components.'^'^ In 
addition, entropic factors may contribute to the monomer incompatibility, which 
are, e.^., related to shape or stiffness disparities.'^^^SQl Since the enthalpic and 
entropic contributions evolve differently as a function of the temperature, the 
X-parameter will in general exhibit a complicated temperature dependence. 




z-2 



(3.15) 



Vig.a{NB = Na)^ 
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Equation-of-state effect : In general, the volume per polymer depends on the com- 
position of the blend. Already the volume per monomer is usually different for 
different monomer species; at constant pressure, x therefore varies roughly lin- 
early with ^A-"^^ In blends of monomers with very similar monomer structure, 
e.g., isotopic blends, the linear contribution vanishes and a weak parabolic de- 
pendence remains, which can partly^^ (but not fulljS^l^ \)q explained by an 'excess 
volume of mixing'. 

Chain correlations : Since the demixing is driven by intermolecular contacts, in- 
tramolecular contacts only contribute indirectly to the ^-parameter. The estimate 
for X can be improved if one replaces the factor z — 2 by an 'effective coordina- 
tion number' z^g which is given by the mean number of interchain contacts per 
monomer. Moreover, the ideality assumption (see Sec. 111.4.2^ is not strictly 
valid. Chains in the minority phase {e.g., A chains in a B-rich phase) tend to 
shrink in order to reduce unfavorable contacts.^S Since the segregation is effec- 
tively driven by x^: X slightly depends on the chain length as a result. The 
situation becomes even more complex if copolymers are involved, which assume 
dumbbell shapes even in a disordered environment.'^'^ This also may affect the 
effective x-parameter.'^Sl 

Composition correlations : The effective interactions between monomers change if 
the local environment is not random. We have already noted earlier that the 
composition may fluctuate. Large-scale fluctuations can be incorporated in the 
Flory-Huggins framework in terms of a fluctuating field theory (see Sees. IIII.2.5I 
and lIII.3.ip . Fluctuat ions (correlations) on the monomeric scale renormalize the 
X-parameter {nonrandom mixing). Moreover, the local fluid structure may depend 
on the local composition {nonrandom packing). 

In view of these complications, establishing an exhaustive theory of the x-parameter 
remains a formidable taskW^^ Even the reverse problem of designing simplified 
particle-based polymer models with a well-defined x-parameters turns out to be highly 
non- trivial.'53 The very concept of a x-parameter has been challenged repeatedly, e.g., 
by Tambasco et a/.,'^ who analyzed experimental data for a series of blends and found 
that their thermodynamic behavior can be related to a single '17 — 1-parameter', which is 
independent of composition, temperature and pressure. They suggest that this param- 
eter may be more appropriate to characterize blends than the x-parameter. However, it 
has to be used in conjunction with an integral equation theory, the BGY lattice theory 
by Lipson,'S3E2] which is much more involved than the Flory-Huggins theory especially 
when applied to inhomogeneous systems. 

Freed and coworkerd^Sl'^B'SSllSQ] j^g^yg proposed a generalized Flory-Huggins theory, 
the 'lattice cluster theory', which provides a consistent microscopic theory for macro- 
scopic thermodynamic behavior. In a certain limit (high pressure, high molecular 
weight, fully flexible chain), this theory reproduces a Flory-Huggins type free energy 
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with an effective x-paramete: 



Xe. = J, + Xo(l ) , (3.16) 

where xo is given by Eq. ()3.15p and the parameters r-,-, pj depend on the structure 
of the monomers j. Dudowicz et al. have recently pointed out that this model can 
account for a wide range of experimentally observed miscibility behavior.l^l^ 

Another pragmatic and reasonably successful approach consists in using x as a 
heuristic parameter. Assuming that it is at least independent of the chain length, 
it can be determined experimentally, e.g., from fitting small-angle scattering data to 
theoretically predicted structure factors-l^^^^ Alternatively, x can be estimated from 
atomistic simulations.'^'^ The results from experiment and simulation tend to compare 
favorably even for systems of complex polymers.'^ 



III. 2 Self-consistent Field Theory 

I have taken some care to discuss the Flory-Huggins theory because it establishes a 
framework for more general theories of polymer blends. In particular, it provides the 
concept of the Flory-Huggins parameter x-, which will be taken for granted from now 
on, and taken to be independent of the composition, despite the question marks raised 
in the previous section (Sec. IIII.1.3p . In this section, I present a more sophisticated 
mean-field approach, the self-consistent field (SCF) theory. It was first proposed by 
Helfand and coworkers^^ ' '''^ and has since evolved to be one of the most powerful tools 
in polymer theory. Reviews on the SCF approach can be found, e.g., in Refs.'^'^'^ 



III. 2.1 How^ it works in principle 

For simplicity, I will first present the SCF formalism for binary blends, and discuss 
possible extensions later. Our starting point is the Edwards Hamiltonian for Gaussian 
chains, Eq. (13. 4p . with a Flory-Huggins interaction term, 

'Hi[pA,PB]/kgT = POX j dr^A^B (3.17) 

where we have defined the 'monomer volume fractions' = Pj/po and incompressibility 
is requested, i.e., ^a + = 1 everywhere. The quantity pj depends on the paths 
of chains of type j and has been defined in Eq. (j3.10p . 

We consider a mixture of ua homopolymers A of length Na and ub homopolymers 
B of length Nb in the volume V. The canonical partition function is given by 

e-Poxfdr^A^s (3.18) 

The product a runs over all chains in the system, pqJ DRq denotes the path inte- 
gral over all paths Ro(s), and we have introduced the factor po in order to make Z 
dimensionless. 



Z 



uaIubI 
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The path integrals can be decoupled by inserting delta- functions / Vpj 6{pj — pj) 
(with j =A,B), and using the Fourier representations of the delta-functions S{pj — pj) = 

/ioo^^i e^^^(^^-^^) and 6{pA + pB - po) = Loo'^C ^'■''^+^^-'"'1 Here N is some 
reference chain length, and the factor 1/N is introduced for convenience. This allows 
one to rewrite the partition function in the following form: 

^oc / VWaVWbV^ I VpaVpbc-^/^b'^ (3.19) 

J ioo J oo 



with 

kgT N 



XN dr<^A^B- dr Wj^j (3.20) 

j=A,B 

[dvi{^A + ^B-l)-Y. ^^-^ In(-OO^)}, 



j=A,B 

{^j = pj/po), where Vj denotes the partial volume occupied by all polymers of type j 
in the system, and the functional 

/I N ■ 
-p^^-nam/k^T g-^/„^d.w,(R(s))^ (3^21) 

is the partition function of a single, noninteracting chain j in the external field Wj. 
Thus the path integrals are decoupled as intended, and the coupling is transferred to 
the integral over fluctuating fields Wj and ^. 

Now the self-consistent field approximation consists in replacing the integral (I3.19P 
by its saddle point, i.e., minimizing the effective Hamiltonian Ti with respect to the 
variables Pj{r) and Wj{r). The minimization procedure results in a set of equations, 

i^j) = with j = A,B 

Wj = x^i^i)-^ with i,j = A,B and i^j, ^' ' 

where {^j) denotes the average of ^j in a system of noninteracting chains subject 
to the external fields Wj{r). One should note that the latter are real, according to 
Eq. (|3.22p . even though the original integral (j3.19|) is carried out over the imaginary 
axis. Intuitively, the Wj{r) can be interpreted as the effective mean fields acting on 
monomers due to the interactions with the surrounding monomers. Together with the 
incompressibility constraint, the equations (j3.22p form a closed cycle which can be 
solved self-consistently. 

For future reference, we note that it is sometimes convenient to carry out the saddle 
point integral only with respect to the variables Wj{r). This defines a free energy 
functional J-scf [^a] ; which has essentially the same form as J-' (Eq. I3.20p , except that 
the variables Wa(r) and ^(r) are now real Lagrange parameter fields that enforce {^a) = 
^A, {^b) = ^b, and $b = 1 — ^A and depend on The same functional can also 
be derived by standard density functional approaches, using as the reference system a 
gas of non-interacting Gaussian chains.!^ 
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In some cases, one would prefer to operate in the grand canonical ensemble, i.e., 
at variable polymer numbers nj. The resulting SCF theory is very similar. The last 
term in Eq. ()3.20p is replaced by (— ZjQj)-, where Zj is proportional to the fugacity 
of the polymers j. 

The formalism can easily be generalized to other inhomogeneous polymer systems. 
The application of the theory to multicomponent A/B/C/... homopolymer blends 
with a more general interaction Hamiltonian "HiIpa, Pb, Per "] replacing Eq. (j3.17p is 
straightforward. The self-consistent field equations (j3.22p are simply replaced by 

If copolymers are involved, the single-chain partition function Qc for the correspond- 
ing molecules must be adjusted accordingly. For example, the single-chain partition 
function for an A:B diblock copolymer of length Nc with A-fraction / reads 

The general SCF free energy functional for incompressible multicomponent systems is 
given by 

-^scf[{p,}] ni[{p,}] If 1 

-J]n, ln(/,o— ), (3.25) 

where the sum j runs over monomer species and the sum a over polymer types. 

The SCF theory has been extended in various ways to treat more complex sys- 
tems, e.g., compressible melts and solutions macromolecules with complex archi- 
tecture,'^ semiflexible polymer^^ with orientational interactions ,'2S1ES1 charged poly- 
merSjl^ polydisperse systems,'^'^ polymer systems subject to stresses,'^'^ systems of 
polymers undergoing reversible bonds,'S51ISSl or polymer /colloid composites.lSStSSl 

III. 2. 2 How it works in practice 

In the previous section, we have derived the basic equations of the SCF theory. Now 
I describe how to solve them in practice. The first task is to evaluate the single-chain 
partition functions Qj and the corresponding density averages {pj) for noninteracting 
chains in an external field. 

We consider a single ideal chain of length N in an external field W{r, s), which may 
not only vary in space r, but also depend on the monomer position s in the chain in 
the case of copolymers. It is convenient to introduce partial partition functions 

q{r, s) = jvR e-^GMA^T !^ ds' W{n{s'),s') ^(j^^^^ _ (g ^g) 

q\r, s) = jvR e-^GMAs^ Jo ds' W{n{s'),N-s') ^(j^^^^ _ (3 27) 
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where the path integrals PR are carried out over paths of length s. According to the 
Feynman-Kac formulaj^ the functions q{r,s) and q'^{r,s) satisfy a diffusion equation 

-9(r,s) = (-A-^^)q{r,s) (3.28) 
Q-^qHr,s) = ij^-^^ -)qHr,s) (3.29) 

with the initial condition g(r, 0) = g^^(r,0) = 1. Numerical methods to solve diffusion 
equations are available, hence q and q^ are accessible quantities. The single-chain 
partition function can then be calculated via 

Q = jdrq{r,N) = jdrq'<{r,N) (3.30) 

and the distribution of the sth monomer in space is g(r, s) g^(r, — s)/Q. 

More specifically, to study binary homopolymer blends, one must solve the diffusion 
equations for the partial partition functions qj = in an external field W = Wj (with 
j = A,B). The averaged volume fractions of monomers j are then given by 

(cl>,(r)) = - -f / dsqj{r,s)q]{r,N -s). (3.31) 

Po Jo 

in the canonical ensemble, and 

(l>,(r)) = ^ dsq,{r,s)q]{r,N-s) (3.32) 

in the grand canonical ensemble. If AB diblock copolymers with fraction / of A- 
monomers are present, one must calculate the partial partition functions qc and ql in 
the external field W{r, s) = Wa for s < fN^ and W{r) = Wb for s > fN^. The 
contribution of the copolymers to the volume fraction ($a) is nc/poQc /A(r) with 
/yi(r) = Jq^" ds gc(r, s)(?^(r, — s) in the canonical ensemble, and Zc/NfA{r) in the 
grandcanonical ensemble. 

With this recipe at hand, one can calculate the different terms in Eqs. (j3.22p . 
The next problem is to solve these equations simultaneously, taking account of the 
incompressibility constraint. This is usually done iteratively. I refer the reader to 
Section 3.4. in Ref . [9T] for a discussion of different iteration methods. 



III. 2. 3 Application: Diblock Copolymer Blends, Part I 

To illustrate the power of the SCF approach, I cite one of its most spectacular suc- 
cesses: The reproduction of arbitrarily complex copolymer mesophases. In a series 
of seminal papers, Matsen and coworkers have calculated phase diagrams for diblock 
copolymer melts.'^'^ Fig. 13.51 compares an experimental phase diagram due to Bates 
and coworkerd^^^ with the SCF phase diagram of Matsen and Bates.'^ The SCF the- 
ory reproduces the experimentally observed structures. At high values of ('strong 
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Figure 3.5: Experimental phase diagram for polystyrene-polyisoprene diblock copoly- 
mer melts (left, after Refs. [93Vl96p compared with phase diagram obtained with the 
SCF theory by Matsen and coworkers (right, after Ref. [93|) in the coordinates ^-nd 
A block volume fraction /. The labels S,C,L, and G correspond to the structures shown 
in Fig. 13. 3[ In addition, the SCF phase diagram features a close-packed sphere phase 
Sep. All phase transitions are first order except for the disordered/lamellar-transition 
at / = 1/2 in the SCF phase diagram. Courtesy of Mark W. Matsen, adapted from 
M.W. Matsen, J. Phys.: Cond. Matter 14, R21 (2002). 



segregation') the SCF phase diagram features the correct sequence of mesophases at 
almost the correct value of the fraction of A- monomers /. At low values of ('weak 
segregation'), the two phase diagrams are distinctly different. This can be explained by 
the effect of fluctuations and will be discussed further below (Sees. [111.2. 51 and IIV.4.21 
see also Fig. ES])- 

Tyler and Morse have recently reconsidered the SCF phase diagram and predicted 
the existence of yet another mesophase, which has an orthorombic unit cell and an Fddd 
structure and intrudes in a narrow regime at the low x-^-end of the gyroid phase.'^ 
This phase was later indeed found in a polystyrene-polyisoprene diblock copolymer 
melt by Takenaka and coworkers .'SS 



III. 2.4 Related Mean-Field Approaches 

So far, I have focussed on sketching a variant of the SCF theory which was originally 
developed by Helfand and coworkers.!^ A number of similar approaches have been 
proposed in the literature. 

Scheutjens and Fleer have developed a SCF theory for lattice models,'^ which is 
applied very widely Scheutj ens-Fleer calculations are very efficient and incorporate in 
a natural way the finite (nonzero) range of monomer interactions. To account for this 
in the Helfand theory, one must introduce additional terms in ()3.9p .^^ which indeed 
turn out to become important in the vicinity of surfaces.*^ 

Carignano and Szleifei'^^have proposed a SCF theory where chains are sampled as a 
whole in a surrounding mean field. Hence intramolecular interactions are accounted for 
exactly and the chain statistics corresponds to that of self-avoiding walks (Sec. III.4.1[) . 
This approach is more suitable than the standard SCF theory to study polymers in 
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solution, or melts of molecules with low-molecular weight, where the ideality assumption 
(see Sec. III.4.2P becomes questionable.'^^ 

In this chapter, I have chosen a field-theoretic way to present the SCF theory. Freed 
and coworker^^l^^ have derived the same type of theory from a density functional 
approach, using a reference system of non-interacting Gaussian chains. Compared to 
the density functional approach, the field-theoretic approach has the advantage that 
the effect of fluctuations can be treated in a more transparent way (see Sec. 1111.2.5^ . 
On the other hand, information on the local liquid structure of the melt, (i.e., monomer 
correlation functions, packing effects etc.), can be incorporated more easily in density 
functional approaches .HQSlUQll Density functionals have also served as a starting point for 
the development of dynamical theories which allow to study the evolution of multiphase 
polymer blends in time^^-^^ (see Sec. IIV.3.2p . 

III.2.5 Fluctuation Effects 

Mean-field approaches for polymer systems like the SCF theory tend to be quite success- 
ful, because polymers overlap strongly and have many interaction partners. However, 
there are several instances where composition fluctuations become important and may 
affect the phase behavior qualitatively. 

To illustrate some of them, we show the phase diagram of ternary mixtures con- 
taining A and B homopolymers and AB diblock copolymers in Fig. 13.61 The left graph 
shows the experimental phase diagram ,1^^ the right graph theoretical phase diagrams 
obtained by Diichs et a|il2Jliiii from the SCF theory (solid lines) and from field-based 
computer simulations (dashed line, see Sec. lIV.2^ for details on the simulation method). 
Regions where different types of fluctuations come into play are marked by I - IV. 

I) Fluctuations are important in the close vicinity of critical points, i.e., continuous 
phase transitions. They affect the values of the critical exponents, which charac- 
terize e.g., the behavior of the specific heat at the transition.!^ In Fig. 13. 6( such 
critical transitions are encountered at high homopolymer concentration, where 
the system essentially behaves like a binary A/B mixture with a critical demix- 
ing point. This point belongs to the Ising universality class, hence the system 
should exhibit Ising critical behavior. It has to be noted that in polymer blends, 

critical exponents typically remain mean-field like until very close to the critical 
point.ESHnil 

II) The effect of fluctuations is more dramatic in the vicinity of order-disorder tran- 
sitions (ODT), e.g., the transition between the disordered phase and the lamellar 
phase at low homopolymer concentrations. Fluctuations destroy the long-range 
order in weakly segregated periodic structures, they shift the ODT and change 
the order of the transition from continuous to first order (Brazovskii mecha- 
nismfnSHnSl) . This effect accounts for the differences between the experimental 
and the SCF phase diagram in Fig. 13.51 



III) 



The SCF phase diagram features a three-phase (Lamellar + A + B) coexistence 
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Figure 3.6: Experimental phase diagram for symmetric ternary blends of 
PDMS+PEE+PDMS-PEE {Nc ~ 5Nh) featuring lamellar phase, phase-separated re- 
gion 2<I>, and microemulsion channel BfiE (left, from Ref . 11091 ). compared with theoret- 
ical phase diagrams (right) from SCF theory (solid lines) and Monte Carlo simulations 
at C = 50 (dashed lines/circles). The dotted lines separate disordered regions with dif- 
ferent local structure: D^E, defect driven disorder, and Gfj,E, genuine microemulsion 
morphology. The regions I-IV are discussed in the text. Left Figure: Reprinted with 
permission from M.A. Hillmyer, W.W. Maurer, T.P. Lodge, F.S. Bates, K. Almdal, J. 
Phys. Chemistry B 103, 4814 (1999). Reproduced by permission of the Royal Society 
of Chemistry. Right Figure: Adapted from Refs. dlOKIII] (Diichs et al 2003 and 2004). 



region reaching up to a Lifshitz point. Lifshitz points are generally believed to 
be destroyed by fluctuations in three dimensions. 

IV) In strongly segregated mixtures, fluctuations affect the large-scale structure of 
interfaces. Whereas mean-field interfaces are flat, real interfaces undulate. The 
so-called 'capillary waves' may destroy the orientational order in highly swollen 
lamellar phases. A locally segregated, but globally disordered 'microemulsion' 
state intrudes between the homopolymer-poor lamellar phase and the homopolymer- 
rich two-phase region in Fig. 13.61 

Both in the cases of III) and IV) , the effect of fluctuations is to destroy lamellar order 
in favor of a disordered state. However, the mechanisms are different. This is found to 
leave a signature in the structure of the disordered phase, which is still locally structured 
with a characteristic wavevector g*.liiiJ In the Brazovskii regime, the wavevectore q* 
corresponds to that calculated from the SCF theory (defect driven disorder regime, 
DjiE). In the capillary wave regime, the characteristic length scale increases, compared 
to that calculated from the SCF theory (genuine microemulsion regime, Gfj,E). 

Formally, the effect of fluctuations is hidden in the overall prefactor po/N in the 
SCF Hamiltonian T-L (Eq. I3.25p . The larger this factor, the more accurate is the sad- 
dle point integration that lies at the heart of the SCF approximation. One can thus 
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define a 'Ginzburg parameter' C = RgPo/N, which characterizes the strength of the 
fluctuations. Here the factor Rg must be introduced to make C dimensionless {d is 
the spatial dimension), and Rg = Nb'^ /6 is the natural length scale of the system, the 
radius of gyration of an ideal chain of length N. The Ginzburg parameter roughly 
corresponds to the ratio of the volume spanned by a chain, Rg, and the volume ac- 
tually occupied by a chain, N/po, and thus measures the degree of interdigitation of 
chains. At C — )• oo, the SCF approximation becomes exact. The numerical simulations 
shown in Fig. 13.61 were carried out at C = 50 (using the length of the copolymers as 
the reference length), which still seems large. Nevertheless, the effect of fluctuations is 
already quite dramatic. 

In three dimensions (d = 3), the Ginzburg parameter is proportional to the square 
root of the chain length, \/iV. This is why the mean-field theory becomes very good 
for systems of polymers with high molecular weight, and only fails at selected points in 
the phase diagram. The relation C oc ^/N has motivated the definition of an 'invariant 
polymerization index' N = Nb^pQ oc C^, which is also often used to quantify fluctuation 
effects.'^'^ In two dimensional systems, C is independent of the chain length and 
fluctuation effects are much stronger. Furthermore, topological constraints become 
important {i.e., the fact that chains cannot cross each other) which are not included in 
the Helfand model. In three dimensional systems of linear polymers, they only affect 
the dynamics (leading to reptation), but in two dimensions, they also change the static 
properties qualitatively.^^ 

Finally, we note that fluctuations can also be treated to some extent within the 
SCF theory, by looking at Gaussian fluctuations about the SCF solution.lH^ This is 
useful for calculating structure factors and carrying out stability analyses. However, 
Gaussian fluctuations alone cannot bring about the qualitative changes in the phase 
behavior and the critical exponents which have been described above. 

III. 3 Analytical Theories 

The SCF equations have to be solved numerically, which can be quite challenging from 
a computational point of view. In addition, they also serve as a starting point for the 
derivation of simpler approximate theories, which may even have analytical solutions 
in certain limits. 

Two main regimes have to be distinguished here. In the weak segregation limit, 
xN is small, and the A and B homopolymers or copolymer blocks are barely demixed. 
This is the realm of the 'random phase approximation' (RPA), which can be derived 
systematically from the SCF theory. In the strong segregation limit, is large, the 
polymers or copolymer blocks are strongly demixed and the system can basically be 
characterized in terms of its internal interfaces. 

III. 3.1 Weak Segregation and Random Phase Approximation 

We first consider the situation at low x^- In this case, the composition varies smoothly, 
A-rich domains still contain sizeable fractions of B-monomers and vice versa, and the 
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interfaces between domains are broad, i.e., tlieir width is comparable to the radius of 
gyration of the chains. 

The idea of the RPA is to perform a systematic expansion about a homogeneous 
reference state. More precisely, we use the SCF free energy density functional, Eq. 
(|3.25p . as a starting point, and then expand J- about the homogeneous state. Defining 

= Pa/po as usual, using + ^ ^, and introducing the Fourier representation 
<I'(k) = J dr e*''''<5(r), we obtain a functional of the form 

^RPA[^A(r)] = y/,oAo„.o + A:,r^{^ j;|c^^(k)|2r2(k) 

^ ^A(k)$A(k')^A(-k-k')r3(k,k') + •••}, (3.33) 

k,kYO 

where /homo is the SCF free energy per chain in the reference system, and the coefficient 
r„ depend on the direct monomer interactions and on the intrachain correlations of 
free ideal Gaussian chains. 

We focus on the leading coefficient. To calculate r2 for a given blend, we define the 
pair correlators i^j^lHJlliil 

K,,{k) = {p^{k)p,{k)) (3.34) 

which give the density-density correlations in an identical blend of noninteracting, ideal 
Gaussian chains and can thus be expressed in terms of Debye functions gjyx) (Eq. (|3.7p ). 
For example, for binary blends of homopolymers with chain length Nj, gyration radii 
Rg^j, and mean volume fractions (j=A,B), the pair correlators are given by 



KAA = —^Ago{k''RlA), KBB = ^^B9oif^^RlB), Kab = KbA = (3.35) 



^^A9oik'RlA), KBB = f^ 

For pure diblock copolymer blends with A-fraction /, one gets 



Kaa = fg^if k^Rl) Kbb = (1 - f?9^{{l - f) k^Rl) 
Kab = k{9oik'Rl) - Pguifk'RD - (1 - f?9oii^ " f)k'Rl 

Having calculated Kij, one can evaluate r2 according tcpl"^^ 

/ Kaa + Kbb + Kab + Kba 



(3.36) 



2xN]. (3.37) 



^ KaaKbb - KabKba 

The function r2 is particularly interesting, because it is directly related to the structure 
factor of the homogeneous phase, ^(k) oc r2(k)~^. Hence the RPA provides expressions 
for structure factors which can be compared to small angle scattering experiments, 
e.g., to determine effective interaction parameters. This is probably its most important 
application. 

We will now discuss specifically the application of the RPA to binary homopolymer 
blends and to diblock copolymer blends. 
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(i) Binary homopolymer blends and Flory-Huggins-de Gennes functional 

According to Eqs. (|3.35p and (|3.37p . the RPA coefficient r2 for binary blends is 
given by 

We assume that the composition varies only slowly in the system (on length scales not 
much shorter than Rg), and expand r2 for small wave vectors. Using ^^(x) ~ 1 — x/3 
and R'^ j = b'jNj/6, and inserting our result in the RPA expansion (j3.33p . we obtain 
the free energy functional 



k2 u2 



+ — y -(^ + ^)|<I>A(k)|2 k"]. (3.39) 

kT^O 

The first two terms in (|3.39p correspond to the second order expansion of the integral 
po / dr /scF(^yl)5 where ^scvi,'^ A) is the SCF free energy per chain in a homogeneous 
system with A- volume fraction ^a- It thus seems reasonable to replace them by the full 
integral. The last term is a square gradient term in real space. Together, one recovers 
the Flory-Huggins-de Gennes functional of Sec. IIII.Ol Eq. (|3.13p . 

(ii) Copolymer melts, Leibler theory, and Ohta-Kawasaki functional 

In diblock copolymer blends, Eqs. (|3.36p and (|3.37p yield the RPA coefficient 

^'^^^ = (g(/)G(1 - /) - (G(l) - G(/) - G(l - /))2/4 " ^^^) ^^-^^^ 

with the short hand notation G(/) = /^^^(//c^iig). At low xA^, r2(k) is positive. 
Upon increasing ■, oiis encounters a spinodal line where r2(k) becomes zero for 
some nonzero k = q*, and the disordered state becomes unstable with respect to an 
ordered microphase separated state. 

Since the function r2(k) is spherically symmetric in k, it does not favor a specific 
type of order. The information on possible ordered states is contained in the higher 
order coefficients r„, most notably, in the structure of the cubic term, T^. In a seminal 
paper of 1979, Leibler has carried out a fourth order RPA expansion and deduced a 
pha se d iagram which already included the three copolymer phases L, C, and S (Fig. 
[3j3j)fnHl Milner and Olmsted later showed that the Leibler theory is also capable of 
reproducing the gyroid phase.'I^ The RPA phase diagram roughly coincides with the 
full SCF phase diagram, as established 1996 by Matsen and Bates,!^ up to < 
12. Unfortunately, fluctuations have a massive effect on the phase diagram at these 
small x^ (see Fig. 13. Sp . therefore the predictive power of the Leibler theory must 
be questioned. Nevertheless, it is useful for identifying potential ordered phases and 
phase transitions in copolymer systems. Generalized Leibler theories still prove to be 
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efficient tools to analyze phase transitions in complex copolymer blends by analytical 
considerations 

Next we attempt to construct a simplified free energy functional for diblock copoly- 
mer melts, in the spirit of the Flory-Huggins-de Gennes functional. To this end, we 
again expand r2 in powers of k, as in (i). Compared to homopolymer blends, how- 
ever, there is an important difference: r2(k) has a singularity at k — )■ and diverges 
according to 

The singularity accounts for the fact that large-scale composition fluctuations are not 
possible in copolymer blends, since the A- and B- blocks are permanently linked to 
each other. It ensures that the structure factor S'(k), vanishes at k — t- 0, suppresses 
macrophase separation and is thus ultimately responsible for the onset of microphase 
separation in the RPA theory. 

A 1/P term like ^AT\\ in a density functional corresponds to a long-range Coulomb 
type interaction. This observation motivated Ohta and KawasakP^D in 1986 to pro- 
pose a free energy functional for copolymer melts, which combines a regular square- 
gradient functional accounting for direct short-range interactions with a long-range 
Coulomb term accounting for the connectivity of the copolymers. In real space, the 
Ohta-Kawasaki functional has the form 



OK 



[^a] = ^ y"dr{w(<I>A) + |(VcI>^)'} + ^^ y"drdr'g(r,r')5$A(r)(5$A(r') (3.42) 



with 5^A{r) = <I>A(r) — ^A- The last term introduces the long-range interactions, with 
G(r, r') defined such that 

Ag(r,r') = -5(r-r'), (3.43) 

which corresponds to ^(r, r') ~ l/|r — r'| in infinitely extended systems. 

Given Eq. (fOT]) . it seems natural to identify A = 3/(2/^(1 - ffR'^g). The choice 
of W and B is somewhat more arbitrary. The function W($yi) is a free energy density 
with two degenerate minima and can be approximated by a fourth order polynomial in 
As for the coefficient of the square gradient term, B, Ohta and Kawasaki originally 
estimated it from the asymptotic behavior of r2 at /c — )• oo, 

r2(k)« 2/(r=7)^'^^ at A;^oo, (3.44) 

which yields B = i?^/(2/(l — /)). Later, they noted that this choice of B gives the 
wrong interfacial width at stronger segregation, which has implications for the elastic 
constants and the equilibrium period of the ordered phases, and suggested to replace 
-B by a constant in the strong segregation limit .'^^D 

The Ohta-Kawasaki functional reproduces microphase separation and complex copoly- 
mer phases such as the gyroid phas^^^HIEl] ^^^^ even the Fddd phase.12il It can be han- 
dled much more easily than the Leibler theory or the full SCF theory (see Sec. IIV.3^ . 
therefore it is particularly popular in large-scale dynamical simulations of copolymer 
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melts (see Sec. IIV.3.2p . Different authors have generalized it to ternary blends con- 
taining cop olymer s .I1231[125]ii26j particular, Uneyama and Doi have recently proposed 
a general density functional for polymer/copolymer blends that reduces to the Flory- 
Huggins-de Gennes functional in the homopolymer case and to the Ohta-Kawasaki 
functional in the diblock case. 

III. 3. 2 Strong Segregation 

I turn to discussing the situation at high x^- The A-rich and B-rich (micro)phases 
are then well-separated by sharp interfaces. The free energy contribution from the 
interfacial regions (i) and the chain conformations inside the A- or B- domains (ii) can 
be treated separately. 

(i) Interfacial profiles and ground state dominance 

In the interfacial region, the free energy is dominated by the contribution of the 
direct A-B interactions and the local stretching of segments. Chain end effects can be 
neglected. This simplifies the situation considerably. 

We first note that the diffusion equation ( (|3.28p or (j3.29p ) for a j-chain or a j-block 
of a chain has the same structure than the time-dependent Schrodinger equation, if 
one identifies s -f-T- it. As is well known from quantum mechanics, the general solution 
can formally be expressed qj{r,s) = '^«V'n,j(r)e~'^"'^^, where {ipnjiT^), ^nj} 

are Eigenfunctions and Eigenvalues of the operator (6|/6 A — Wj{r)/N). At large s, 
the smallest Eigenvalue eoj dominates, i.e., qj{r,s) cx 'ipQj{r)e~'^°'^'^ , and the resulting 
density in the large Nj limit is pj oc iV'OjPe"'^'''^^^ . This type of approximation is called 
'ground state dominance'. It is commonly used to study polymers at interfaces and 
surfaces. 

In the case of blends, we have the freedom to shift the fields Wj{r) by a constant 
value, hence we can set eoj = 0. The self-consistent field equations can thus be written 
as 

p,=Po\^PJirr with (-i-A-^)V,=0 and W, = ^ - 1 (3.45) 

where ^pj is normalized such that Qj = j dr \ipj(r\'^ = Vj is the partial volume occupied 
by the polymers j, and ^(r) ensures Ylj IV'jP = 1- 

In order to derive an epression for the free energy, we first note that Eqs. (j3.45p 
minimize a Lagrange action, 

C = ni + ^Po Jdr (bliViJAf + 6|(VVb)'), (3.46) 

with respect to Tpj under the constraint iV'Ap + I^Bp = 1- One easily checks that 
C vanishes for homogeneous bulk states, and that the minimized C is equal to the 
extremized SCF Hamiltonian J-scf, Eq. ()3.25p up to a constant. Hence C can be 
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identified with the interfacial free energy. Rewriting it in terms of the volume fractions 
4>j and using [V^a)^ = (V<I'_b)^, one obtains the free energy functional 

TU^Air)] = + Po / dr + ^){V^Af, (3.47) 

which reproduces Eq. (j3.14p for Flory-Huggins interactions ()3.17p . 

For bA = bs = b and Flory-Huggins interactions, the self-consistent field equations 
(I3.45P are solved by a tanh profile, {pA — Pb) ~ Po^siah.{z/wssh) with the interfacial 
width tfssL = ^/\/6x and the interfacial tension crggL = k^p^b^J y^jQ^^^^ 

(ii) Copolymer conformations and strong stretching theory 

The free energy functional (13.470 is sufficient to describe strongly segregated ho- 
mopolymer blends. In copolymer blends, additional contributions come into play due to 
the fact that the copolymer junctions are confined to the interfaces and the copolymer 
blocks stretch away from them into their respective A or B domains. The associated 
costs of configurational free energy can be estimated within a second approximation 
scheme, the 'strong stretching' theory (SST).'I2sHl3i| 

The main idea of the SST was put forward in 1985 by Semenov,1^ who noted 
that for strongly stretched copolymer blocks, the paths fluctuate around a set of 'most 
probable paths'. This motivates to approximate the single-chain partition function Q, 
Eq. (j3.2ip . by its saddle point, i.e., the path integral in Q is replaced by an integral 
over 'classical' paths Rc that extremize the integrand and thus satisfy the differential 
equation^^ 

We will treat the copolymer blocks as independent chains of length M . The classical 
paths corresponding to one block are then characterized by their boundary conditions, 
R'c(O) = Yj and Rc(M) = r^, where the junction r^ is confined to an interface and the 
free end is distributed everywhere in its domain. 

Next, we note that for infinitely long blocks, the classical paths must satisfy 

dRc/ds|s=A/ = for M ^ oo (3.49) 

at the free end. Mathematically speaking, they would not have a well-defined end 
position otherwise. Physically speaking, the 'average' chain representing the classical 
path does not sustain tension at the free end, which seems reasonable. In the following, 
Eq. (I3.49P is also imposed for finite (large) blocks as an additional boundary condition. 
Eq. (j3.48p is then overdetermined and can, in general, no longer be solved for arbitrary 
end positions rg. To ensure that chain ends are indeed free to move throughout the 
domain, the field l^(r) must have a special shape. Specifically, near flat interfaces it 
must be parabolic as a function of the distance z from the interface ll^SJUSOj 
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This is one of the main results of the SST. It generally applies to situations where 
strongly stretched polymers are attached to an interface, e.g., strongly segregated 
copolymer blocksj^'^^ or polymer brushes in solvents of arbitrary quality.'I^ The 
SST field must always have the form (j3.50p . and the remaining task is to realize this 
by a suitable choice of the chain end distribution P{re). In the incompressible blend 
case, -P(re) must be chosen such that the density in the domains is constant, po. 

Luckily, we do not have to evaluate P{ve) explicitly to calculate the free energy. The 
SST field has another convenient property: One can show that the stretching energy of 
classical paths of fixed length in a field satisfying Eqs. ()3.48p and ()3.49p is exactly 
equal to the negative field energy. 

Summing over all blocks in a domain, the total stretching energy is thus given by 

where the integral is over the volume of the domain, and z denotes the closest distance 
to an interface. The total energy of the system can be estimated as the sum over the 
stretching energies in the different domains, Eq. (j3.52p . and the interfacial energy, Eq. 
(|3.47p . and then used to evaluate the relative stability of different phases. In the strong 
segregation limit, only the C, L, and S phase are found to be stable,'^ in agreement 
with SCF calculations at high x-^- 

The validity of the strong stretching theory seems to be restricted to very large 
chains .^^23 This is presumably to a large extent due to the requirement (I3.49p . which 
does not necessarily hold for classical paths of finite length. Netz and Schicli^^'^^ 
have shown that an unrestricted 'classical theory', which just builds on the saddle 
point integration of Q and avoids using (j3.49p . gives results that agree better with 
the SCF theory. However, the classical theory has to be solved numerically, and the 
computational advantage over the full SCF theory is not evident. 

The SST has found numerous application^^ and has been extended and improved 
in various respect. It provides an analytical approach to analyzing multicomponent 
polymer blends in a segregation regime where the SCF theory becomes increasingly 
cumbersome, due to the necessity of handling narrow interfaces. 



III. 4 An Application: Interfaces in Binary Blends 

To close the theory section, I discuss the simplest possible examples of an inhomoge- 
neous polymer system: An interface in a symmetrical binary homopolymer blend. This 
system has been studied intensely in experiments.'^^^'^^^S] gy mixing random copoly- 
mers of ethylene and ethyl-ethylene with two different, but very well defined copolymer 
ratios, Carelh et al^^^ were able to tune the Flory-Huggins parameter very finely 
and study interfacial properties in a wide range of between the weak segregation 
limit and the strong segregation limit. Fig. 13.71 compares the results for the interfacial 
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Figure 3.7: Intrinsic width of interfaces between A- and B-phases in binary polyolefin 
blends as a function of x^i compared with the predictions of the weak segregation 
theory (long dashed line) , the strong segregation theory (dot dashed line) , and the full 
SCF theory (solid line). Prom Ref. [T371 

width and compares them with the mean- field prediction for the weak segregation limit, 
the strong segregation limit, and the full numerical result. 

We should note that there is a complication here. As mentioned earlier, fluid-fluid 
interfaces are never flat, they exhibit capillary waves .'^^HUISI ^phis leads to an apparent 
broadening of the interfacial width The apparent width depends on the lateral 

length scale of the observation L and can be calculated according tJEZlEll 

o n kJT , / L \ o o 2k„T / L\ , , 

u;2 =u>g + ^ ln(— ) or = wl + ^^ ln(—], (3.53) 

depending how it is measured. Here wo is the 'intrinsic' width, a the interfacial tension, 
and ao a 'coarse-graining length', which is roughly given by the interfacial width.'^^SllMSI 
Both the quantities L and oq are not very well defined in an actual experiment. Fortu- 
nately, they only enter logarithmically, therefore the result is not very sensitive to their 
values. The theoretical curves shown in Fig. 13.71 include the capillary wave broadening, 
calculated using the interfacial tension from the respective theory. 

Comparing the curves in Fig. 13. 7|, one finds that the weak segregation theory con- 
sistently overestimates the width, and the strong segregation theory consistently un- 
derestimates it. The numerical SCF values interpolate between the two regimes and 
are in excellent quantitative agreement with the experimental data over almost the 
whole range of x^- The SCF theory is also found to perform very well compared to 
computer simulations .H^ESl n reproduces many features of the interfacial structure, 
such as chain end distributions, local segment orientations etc. at a quantitative level, 
if capillary waves are accounted for,'^^ This illustrates the power of the SCF theory 
to describe the local structure of inhomogeneous polymer systems, even if the global 
structure is affected by large-scale composition fiuctuations. 
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IV Simulations of Multiphase Polymer Systems 

Whereas theoretical work on multiphase polymer systems has a long-standing tradi- 
tion, the field of simulations in this area is much younger. This is because polymer 
simulations are computationally very expensive, which has essentially rendered them 
unfeasible until roughly 20 years ago. In this section, I will attempt to give an overview 
over the current state-of-the-art of simulations of inhomogeneous multicomponent poly- 
mer systems. 

IV. 1 Coarse-Graining, Part II 

One of the obvious challenges in multiphase polymer simulations is that polymers are 
such big molecules, which moreover self-organize into even larger supramolecular struc- 
tures. Polymeric materials exhibit structure on a wide range of length scales, from the 
atomic scale up to micrometers. Their specific material properties are to a great extent 
determined by local inhomogeneities and internal interfaces, and depend strongly on 
the interplay between these mesostructures in space and time. In order to understand 
the materials and make useful predictions for new substances, one must analyze their 
properties on all time and length scales of interest. Therefore, multiscale modeling has 
become one of the big topics in computational polymer science. 

The central element of multiscale modeling is coarse-graining. By successively elim- 
inating degrees of freedom (electronic structure, atomic structure, molecular structure, 
etc.), a hierarchy of models is constructed (see Sec. III.2p . For each type of model, 
optimized simulation methods are developed, which allow one to investigate specific 
aspects of the materials. 

Having identified suitable classes of coarse-grained models, one can proceed in two 
different manners: 

Generic Modeling 

This approach has been favored historically, and up to date, the overwhelming 
majority of simulations of multicomponent polymer systems is still based on it (see 
Sec. IIV.4p . Generic models are simple and computationally efficient. They are not 
designed to represent specific materials; rather, they are the simulation counterpart 
of the theoretical models discussed in the previous section. They are suited to study 
generic properties of polymer and copolymer systems, i.e., to identify the behavior 
that can be expected from their stringlike structure, their chemical incompatibility etc. 
Simulations of generic models are also particularly suited to test theories. 

Generic models are used in all areas of materials science, and in most cases, they only 
give qualitative insights into the behavior of a material. This is different for polymers, 
because of their universal properties (see Sec.[Tl]). For example, we have already seen in 
Fig. 13.71 that a generic theoretical model (the Edwards model) quantitatively predicts 
important aspects of the interfacial structure in real polymer melts. 

Nevertheless, the predictive power of generic models is restricted, and relies on 
the knowledge of 'heuristic' parameters such as the x-parameter. Therefore, a second 
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approach is attracting growing interest. 

Systematic Bottom-Up Modeling 

The idea of systematic coarse-graining is to estabhsh a hierarchy of models for 
the same specific material, starting from an ab-initio description, with well-defined 
quantitative links between the different levels. 

Ideally, the goal is to replace many degrees of freedom by a selection of fewer 'effec- 
tive' degrees of freedom. If one is only interested in equilibrium properties, the problem 
is at least well-defined. For each possible coarse-grained configuration, one must evalu- 
ate a partial partition function of the full system under the constraints imposed by the 
values of the coarse-grained degrees of freedom. This procedure results in an effective 
potential in the coarse-grained space, which is, in general, a true multibody potential 
- it cannot be separated into contributions of pair potentials. If one is interested in 
dynamical properties, the situation is even more complicated. One must replace a dy- 
namical system for all variables by a lower dimensional system for a subset of effective 
variables. This can be done approximately, e.g., using Mori Zwanzig projector opera- 
tor techniques. The new dynamical system is inevitably a stochastic process with 
memory, i.e., the future time evolution not only depends on the current state of the 
system, but also on its entire history. 

Obviously, such 'ideal coarse-graining' is not feasible for polymer systems. Instead, 
researchers adopt a heuristic approach, where they first define a coarse-grained model, 
which typically has no memory and only pair potentials, and match the properties of 
the coarse-grained model with those of the fine-grain model as good as they can: The 
model parameters are chosen such that the coarse-grained model reproduces physical 
properties of interest, such as correlation functions or diffusion constants.'^^^^^^^ 

Already early on, researchers have started to develop schemes for mapping real poly- 
mers on lattice models. ^^^"^^^ Nowadays, off-lattice models are more common. Early 
approaches focussed on the task of reproducing the correct intrachain correlations by 
optimizing the bond potentials in the chains.'^^ Later, the interchain correlations were 
considered as well, which can be matched by adjusting the non-bonded, intermolecular 
potentials in the coarse-grained model. It is important to note that the resulting ef- 
fective potentials depend on the concentration and the temperatur^^SU (much like the 
X-parameter itself). Different methods to determine effective potentials have been de- 
vised and even automated packages are 

available.HSaCSM The reverse problem - how 
to reconstruct a fine-scale model from a given coarse-scale configuration - has also been 
addressed .l^^l^^ Nowadays, the available techniques for mapping static properties are 
relatively advanced. In contrast, the field of mapping dynamical properties is still in 
its infancy.'^^ 

The standard multiscale approach is sequential, i.e., numerical simulations are car- 
ried out separately for different levels. Currently, increasing effort is devoted to devel- 
oping hybrid schemes where several coarse-graining levels are considered simultaneously 
within one single simulation.'^^'I^ 

Despite the large amount of work that has already been devoted to systematic 
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coarse-graining, coarse-grained simulation studies of realistic multicomponent poly- 
mer blends are still scarce. Mattice and coworkers have carried out lattice simula- 
tions of blends containing polyethylene and polypropylene homopolymers and copoly- 
j^gj,g|i68lti70i pg^iigj- Q,/ have developed and studied a coarse-grained model for blends 
of polyisoprene and p olystyr ene .HZUSZS] 

IV. 2 Overview of structural models 

After these general remarks, I shall give a brief overview over the different models that 
are currently used in multicomponent polymer simulations. 

IV.2.1 Atomistic Models 

Atomistic simulations are computationally intensive, and rely very much on the quality 
of the force fields. (Force fields are a separate issue in multiscale modeling, which shall 
not be discussed here.) Therefore, atomistic simulations of blends are still relatively 
scarce. So far, most studies have focussed on miscibility aspects.l^^^^^^^ Already early 
on, atomistic and mesoscopic simulations were combined in multiscale studies: Atom- 
istic simulations were used to determine the Flory-Huggins x-parameter, coarse-grained 
methods were then applied to study large-scale aspects of phase separationP^^'''^ or 
mesophase formation.'!^ Only few fully atomistic studies deal with aspects beyond 
miscibility, e.g, the formation of lamellar structures in diblock copolymers, or the 
diffusion of small molecules in blends. 

IV. 2. 2 Coarse-Grained Particle Models 

The coarse-grained models for polymers can be divided into two main classes: Coarse- 
grained particle models operate with descriptions of the polymers that are considerably 
simplified, compared to atomistic models, but still treat them as explicit individual 
objects. Field models describe polymer systems in terms of spatially varying continuous 
fields. I begin with discussing some of the most common particle models. 

Lattice Chain Models 

Lattice models have the oldest tradition among the coarse-grained particle models 
for polymer simulations, and are still very popular. The first molecular simulations of 
multiphase polymer systems - studies of binary homopolymer blends by Sariban and 
Binder in IOSTSsSIIISSI a,nd by Cifra and coworkers in lOSS^^ - were based on lattice 
models. They are particularly suited to be studied with Monte Carlo methods, and 
several smart Monte Carlo algorithms have been designed especially for lattice polymer 
simulations . 

In molecular lattice models, the polymers are represented as strings of monomers 
confined to a lattice. A natural approach consists in placing the 'monomers' on lattice 
sites and linking them by bonds that connect nearest-neighbor sites. For many appli- 
cations, it has proven useful to apply less rigid constraints on the links and allow for 
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bonds of variable length, which may also connect second-nearest neighbor^^^ol qj- stretch 
over even longer distances Moreover, the lattice is usually not entirely filled with 
monomers, but also contains a small fraction of voids. This is because most Monte 
Carlo algorithms for polymers do not work at full filling, and special algorithms have 
to be devised for that case.'222J One particularly popular lattice model is the 'bond- 
fluctuation model', devised in 1988 by Carmesin and Kremer.l^^si jg leased on the 
cubic lattice; monomers do not occupy single sites, but entire cubes in a cubic lattice. 
They are connected by 'fluctuating bonds' of varying length, (2, \/5, \/6, 3 or ^/W lat- 
tice constants). In the bond-fluctuation model, a polymer system behaves like a dense 
polymer melt already at the volume fraction 0.5. Therefore, it can be simulated very 
efficiently. 

Despite their intrinsically anisotropic character, lattice models are able to reproduce 
most known self-assembled mesophases in copolymer melts, even the gyroid phase in 
diblock copolymers.f^H Nowadays, they are used to study such complex systems as 
ABC triblock copolymer melts confined in cylindrical nanotubes,'^^ which feature a 
rich spectrum of novel morphologies, e.g., stacked disks, curved lamellar structures, 
and various types of helices. 

OfF-Lattice Chain Models 

For many years, only lattice simulations were sufficiently efficient that they could 
be used to study polymer blends at a molecular level. With computers becoming more 
and more powerful, off-lattice chain models become increasingly popular. Compared to 
lattice models, they have the advantage that they provide easy access to forces and can 
also be used in Molecular Dynamics or Brownian Dynamics simulations. They do not 
impose restrictions on the size and shape of the simulation box (in lattice models, the 
box dimensions have to be multiple integers of the lattice constant). The structure of 
space is not anisotropic as in lattice models. Whereas the inherent anisotropy of lattice 
models does not seem to cause problems if the lattice model is sufficiently flexible and 
if the chains are sufficiently long, simulations of shorter chains can be hampered by 
lattice artifacts. 

In bead-spring models, polymers are represented by chains of spherically symmetric 
force centers connected by springs. Numerous variants have been proposed, which differ 
in the choice of the spring potentials (the bonded interactions) and the choice of the 
pairwise interactions between the beads (the non-bonded potentials). 

The simplest choice of spring potential is a harmonic potential. In order to prevent 
chains from crossing each other in dynamical simulations, an anharmonic cutoff on the 
spring length is often imposed. A popular choice is the 'Finitely Extensible Nonlinear 
Elastic' (FENE) potential 

l-PENE(6) = ^dMn(l-^^^), (3.54) 

which reduces to a harmonic spring potential with equilibrium spring length bo at 
b ^ bo, and diverges at |6 — 6o| — t- d. In some applications, the springs are constrained 
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to have fixed lengths - however, this requires the use of special algorithms and changes 
slightly the distribution of bond angles.'^^ In addition, some bead-spring models also 
include bending potentials that allow one to tune the chain stiffness, or even torsional 
potentials. 

The non-bonded interactions drive the segregation of the monomers. As we have 
discussed in Sec. IIII.1.31 both energetic and entropic factors can contribute to making 
monomers incompatible. Many models operate with energetic monomer (in) compatibi- 
lities, but models with entropically driven segregation are also common. As an example, 
we consider one commonly used type of potential, the truncated Lennard-Jones poten- 
tials 

V,^j{r)=e(^{a/r)^^ -{a/rf + C^ for r < (3.55) 

(Vlj = otherwise), where the parameter C is chosen such that V^j(r) is continuous 
at r = rc- If the cutoff parameter rc is larger than 2^^^a (a common choice is r = 
2.5a), the potential has a repulsive core surrounded by an attractive well . In this 
case, energetic incompatibility can be imposed by using species dependent interaction 
parameters Cij with 2eAB < ^aa+^bb- If the cutoff parameter is rc < 2^/^, the potential 
is purely repulsive. In this case, one can still enforce monomer segregation by choosing 
species dependent and non-additive interaction radii aij with 2a ab > ctaa + o'bb- 
The mechanism driving the segregation is the Equation-of-State effect discussed in Sec. 
IIII.1.3[ A simulation model that is based on this idea has been proposed by Grest, 
Lacasse and Kremer in 1996.'^^ It is more efficient than conventional models with 
energetically driven segregation, because the interaction range is shorter. We note that 
high pressures have to be applied to keep the chains together and to drive the demixing. 

In the context of 'Dissipative particle dynamics' (DPD) simulations, (see below) 
it has also become popular to use soft non-bonded potentials without a hard core. A 
typical DPD-potential has the form^'^'' 

Vupo{r) = - {^1- r/rcj for r < rc (3.56) 

(^DPD = otherwise). Demixing is induced by using species dependent parameters 
Vij > with 2vAB > {i^AA + vbb)- The mechanism driving the segregation is again 
related to the Equation-of-State effect - like particles overlap more strongly than unlike 
particles. The DPD simulation method was originally proposed in the context of fluid 
simulations, where every DPD particle supposedly represents a lump of true particles. 
This motivates the absence of hard, impenetrable cores and even a roughly linear 
shape.'^SSl Early on, DPD potentials were also used in polymer simulations. '223 As long 
as topological constraints are not important, the monomer potentials do not need to 
have a hard core (see also the discussion at the end of the next section). 

'Edwards' Models 

A special class of chain models are the Edwards models, which shall be treated 
separately. The idea is to implement directly an Edwards-type interaction %! (Eq. 
()3.9p ) in a particle-based simulation. The molecules are modeled as off- lattice chains 
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as before, and the non-bonded interactions are given by a potential V = Tiiiipi}] that 
depends on local monomer densities pi{r). To complete the definition of the model, one 
must prescribe how to evaluate the local densities. This is most easily done by simply 
counting the monomers in each cell of a grid. Other prescriptions that do not impose 
a grid are also conceivable. 

Edwards models are not yet very common in simulations of multiphase polymer 
systems, but they will very likely gain importance in the future. The basic idea was put 
forward by Zuckermann and coworkers^''^ ^^*^ in 1994 in the context of polymer brush 
simulations. It was first applied to studying microphase separation in block copolymers 
by Besold et al.p^ who also showed that the model produces correct single-chain 
behavior in solutionPI2] (^.g., the chains behave like regular self-avoiding chains). In the 
following, the power of the approach has been demonstrated in a series of impressive 
work by Miiller and coworkerd^^ (see also Sec. llV.S.lh . 

Two notes are in order here. First, it is difficult to impose strict incompressibility 
in particle-based simulations. Instead, dense melt simulations usually operate at finite 
compressibility: One introduces an additional ternJ^ 

^/,comp. = "^{PA + PB - paf (3.57) 

with a high modulus k in the interaction Hamiltonian. Second, chains may overlap 
in the Edwards models. They share this 'problem' with the DPD models introduced 
in the previous section. In three dimensions, the lack of hard core interactions has no 
effect on the static properties of linear p olymer s .I21211213J Topological constraints become 
important in two dimensions,'^ or in melts of cyclic polymers,'2SI or, most notably, 
when looking at dynamical properties. In Edwards models, chains do not entangle, 
and reptation dynamics has to be put in 'by hand'.'^El 

Ellipsoid Model 

In 1998, Murat and Kremer proposed a model that allows to study weakly segre- 
gated polymer blends on the scale of the gyration radius Rg and beyond.l^lS! if details 
of the conformations are not of interest, the polymers can be replaced by single, soft 
particles with ellipsoidal shape. The idea was then pursued mainly by Eurich, Maass, 
and coworkers plIlSISI -^^j^q also suggested to extend the model to diblock copolymers 
by modeling them as dimers to account for their dumbbell shape.l^^ Sambrisky and 
Guenza have recently worked out a microscopic foundation for coarse-graining diblock 
copolymers into such dumbbells, which is based on liquid-state integral equations .12221 

IV.2.3 Field-Based Models 

Field-based models for polymer systems no longer treat the molecules as individual 
objects, but describe them in terms of locally varying fiuctuating fields. Among these, 
the molecular field theories still incorporate some information on the conformations of 
molecules, and the Ginzburg- Landau models focus on the large-scale structure only. 
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Field-Theoretical Models 

Field-theoretical models are the molecular field equivalent of the particle-based 
'Edwards' models discussed above. The idea is to use the starting point of the self- 
consistent field theory, i.e., a field-theoretic expression for the partition function {e.g., 
Eq. (j3.19p ). and to evaluate the functional integrals over fluctuating fields by simulation 
methods. Since some flelds are complex, one is confronted with a sign problem (the 
integrand oscillates between positive and negative). Fredrickson and coworkers have 
demonstrated that the integrals can nevertheless be evaluated in many cases using a 
method borrowed from elementary particle physics, the 'Complex Langevin' simulation 
method.'^'^^'^^ We refer to Ref. M for a general presentation of the method, and Ref. 
12231 for technical details. 

To be more speciflc, let us consider the system discussed in Sec. IIII.2.11 a blend of 
polymers/cop olymers containing two types of monomers A and B, with Flory-Huggins 
interactions (|3.17p . Rather than evaluating the integral over all five fiuctuating flelds 
Wa,B:^: PA,B in Eq. (j3.19p . we reconsider the original partition function, Eq. (j3.18p . 
and decouple the integrals over different paths (polymer conformations) by means of 
a Hubbard-Stratonovich transformation.^" This is possible, because the interaction 
Flory-Huggins interaction '^/[{/Oj}] is quadratic in the densities pi. The result is a 
functional integral over 'only' two fields - one which is conjugate to the total density 
P = PA+ Pb and imaginary, and one which is conjugate to the composition m = PB — PA 
and real: 

Z oc / VW+ I VW- e--^FTs/fcB^ (3 58) 

J ioo J oo 

with 

T^^,[W+,W-] = ^[^ jdrW^- ldrW+-^Vo.^ln{poQJn^)]. (3.59) 

This partition function is the starting point for field-theoretical simulations of incom- 
pressible AB (co)polymer systems. 

The idea to study multiphase polymer systems by direct evaluation of fiuctuating 
field integrals like Eq. (j3.58p was first put forward by Ganesan and Fredrickson in 
2001.12211 They used Complex-Lan gevin simulations to look at fiuctuation effects in 
pure symmetric diblocks. Since then, Fredrickson and coworkers have shown that the 
method can be extended in various ways, e.g., it can deal with external stresse^^lll) ^ 
and it can be applied very naturally to charged p olymers .122611227| Nevertheless, there are 
still problems with it. The theoretical foundations of the Complex Langevin method 
are not fully established. One can show that under certain conditions, it produces the 
correct statistical averages if the system reaches equilibrium.'^ However, the stability 
of the simulations cannot be ensured. Even if the Langevin step size is chosen small, 
a fraction of simulation runs still crashes .l^SSl Such problems are unknown from other 
simulation methods. Since it is generally not widely used and very young in the polymer 
community, the Complex Langevin method still suffers from teething troubles. 

An alternative approach to carrying out field-theoretical simulations in melts has 
been suggested by Diichs et. al^^ (see also Ref. [9T|) . They suggested to treat only 
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the real integral over W- with computer simulation and approximate the problematic 
imaginary integral over W+ by its saddle-point. The integral over W- can be evaluated 
by standard Monte Carlo methods. The results from this partial saddle point approach 
were shown to agree quantitatively with the full Complex Langevin simulation.H^^ Pos- 
sible strategies to improve upon the saddle point integral have been pointed out by 
Baurle and coworker s.l^^SHlSol 

As already pointed out, field-theoretic models are in some sense equivalent to 
particle-based Edwards models (if the latter use discrete Gaussian chain models). In 
field-theoretical simulations, the mean- field limit can be reached very naturally and at 
low computational cost. Hence they are more efficient than particle-based simulations 
at high polymer densities, close to the SCF limit, whereas particle-based simulations 
are superior at lower polymer densities. 

An application of field-theoretic simulations was already shown in Sec. IIII.2.5I The 
simulation data shown in Fig. 13.61 were obtained with field-theoretic Monte Carlo sim- 
ulations. 

Molecular Density Functionals 

Molecular density functionals are free energy functionals of the type (I3.25p . They 
are used in dynamical density functional simulations (see Sec. IIV.3.2p . 

Ginzburg-Landau Models 

Ginzburg-Landau models no longer incorporate specific information on chain confor- 
mations and thus have a much simpler structure than molecular field models. Examples 
are the Flory-Huggins-de Gennes functional for homopolymer blends (Eq. (|3.13p ). or 
the Ohta-Kawasaki functional for diblock copolymer melts (Eq. (j3.42p ). This is the 
highest level of coarse-graining discussed in the present article. 

IV. 3 Overview of Dynamical Models 

In many cases, only static equilibrium properties are of interest, and then most dynam- 
ical models are equivalent. When looking at dynamical properties such as dynamical 
response functions, or at nonequilibrium situations, the choice of the dynamical model 
becomes relevant. In the following, I summarize some important models that have been 
used to study multiphase polymer systems. 

IV. 3.1 Particle-Based Dynamics 

Kinetic Monte Carlo (MC) 

A priori, the Monte Carlo (MC) method has been invented as a method to eval- 
uate high dimensional integrals {i.e., thermal averages), and is designed for studying 
dynamics. Nevertheless, MC simulations are used for dynamical studies, based on the 
fact that like many static properties, dynamical phenomena are also often governed 
by universal principles. In kinetic MC simulations, one analyzes the artificial Monte 
Carlo evolution of a system in order to gain insight into real dynamical processes in 
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the system. The main requirement is that the Monte Carlo moves are only local and 
reasonably 'realistic', i.e., chain crossings are not allowed. Kinetic Monte Carlo simu- 
lations have been used, e.g., to study the early stages of demixing in binary blends or 
the ordering dynamics in block copolymer melts. 

Molecular Dynamics (MD) 

In Molecular Dynamics simulations, one solves directly Newton's or Hamilton's 
equations of motions. This is the most straightforward approach to studying dynamical 
processes in a system. As usual, of course, the devil is in the details.-^^ 

Brownian Dynamics 

In Brownian dynamics simulations, one also solves equations of motion, but the un- 
derlying dynamics is not Hamiltonian. As we have noted earlier (Sec. IIV.T]) . systematic 
coarse-graining necessarily turns a dynamical system into a stochastic process. Brown- 
ian Dynamics simulations account for this fact. They include dissipative and stochastic 
forces, which supposedly represent the interaction of the coarse-grained degrees of free- 
dom with those that have been integrated out. In general, however, these forces are 
not derived systematically from the original model, but postulated heuristically. 

More specifically, the particles j experience three types of forces: 

mjVj = Ff + Ff + Ff (3.60) 

The first term F^* correspondes to the conservative forces, which are derived from 
the interparticle potentials of the structural model under consideration. These forces 
also enter the standard Molecular Dynamics simulations. The second Fj^ term is a 
dissipative force that couples to the velocity of the particles. The last term Fj^ is a 
Gaussian stochastic force with mean zero whose amplitude is related to the dissipative 
force by a 'fluctuation-dissipation theorem'. In a canonical simulation, the last two 
terms constitute a thermostat, i.e., they maintain the system at a given temperature 
T. Microcanonical models where energy is randomly shifted between 'internal' and 
'external' degrees of freedom have been designed as well. 

In the simplest (canonical) Ansatz, the dissipative force on a particle j at time t 
is proportional to its velocity Vj(t). According to the fluctuation-dissipation theorem, 
the stochastic force then fulfills 

Ff = -7,v, ^ (t) F«^(t')) = 27, k,T 5,k 5{t - t'). (3.61) 

Other choices are possible. For example, the dissipative force on j can depend on the ve- 
locities of other particles k as well, and/or on the history of the system. The fluctuation- 
dissipation relation for the stochastic force has to be adjusted accordingly.'^^ The dy- 
namics defined by Eq. (j3.6ip is commonly used because it is so simple, but it has the 
disadvantage that it is not Galilean-invariant - a system moving at constant speed is 
treated differently than a system at rest. Therefore, it does not incorporate hydrody- 
namic effects correctly, and it is not suited to study nonequilibrium systems such as 
polymer melts in shear flow. 
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Dissipative Particle Dynamics (DPD) 

The problems of Eq. ()3.6ip are avoided in the recently developed 'Dissipative Par- 
ticle Dynamics' (DPD) method P^^IISM] j^pj) ig g, special type of Brownian Dynamics 
which is Galilean-invariant and has become very popular in simulations of complex flu- 
ids in recent years. The dissipative and stochastic forces are constructed such that they 
conserve the total momentum and angular momentum in the system. Consequently, 
they couple to relative velocities of particles rather than absolute velocities, and they 
act as central forces (to preserve angular momentum). Specifically, the forces acting 
on a particle j are given by 

Ff - ^ |7tj(rjfc) {rjk^-jk) rjk + k^T uj{rjk) rjk Cjfc}, (3.62) 

where r^fc = rj — is the vector separating two particles j and k, rj^. its length, Tj^ the 
unit vector in the same direction, uj{r) an arbitrary function with a cutoff, and Cjfc are 
symmetric and uncorrelated Gaussian random numbers with mean zero and variance 
one. 

In the literature, the term 'DPD simulations' often refers to simulations that use 
DPD dynamics in combination with soft 'DPD potentials' (see Sec. lIV.2^ (Eq. (I3.56P ). 
However, DPD can of course be used as a dynamical modeP^ or simply as a thermo- 
stat^^^ in combination with any type of potential. 

In contrast to kinetic Monte Carlo simulations or simple Brownian dynamics simu- 
lations (using Eq. (j3.6ip ). DPD simulations take full account of hydrodynamic interac- 
tions. Studies of microphase separation in copolymer melts have shown that this makes 
a difference. The dynamics is strongly affected by hydrodynamic effects in certain re- 
gions of phase space - in particular, hydrodynamic interactions play a critical part in 
helping the system to escape from metastable transient states.'^^Sl 

Single Chain in Mean Field (SCMF) 

In 2005, Miiller proposed an efficient method to study the ordering dynamics of 
polymer blends within Edwards models (see Sec. IIV.2.2p .^'^^^^- The idea is to take 
snapshots of the density configurations in regular intervals, and to let the chains move 
in the fields created by these snapshots, e.g., by kinetic Monte Carlo simulations. If the 
fields were updated after every Monte Carlo move, this would correspond to a regular 
simulation of the Edwards model. Daoulas et al have shown that it is possible to update 
much less oftenP^ without changing the results. This makes the method very efficient 
and especially suited for the use on parallel computers. 

A similar idea had been put forward already in 2003 by Ganesan and Pryamitsin,'^^ 
in a less transparent formulation that involves self-consistent fields, to study stationary 
inhomogeneous polymer systems in an externally imposed fiow field. In the absence of 
fiow, the model of Ganesan and Pryamitsin is equivalent to that of Miiller et al. A hy- 
brid method that combines kinetic Monte Carlo and the self-consistent field formalism 
has also recently been proposed by Baurle and Usami.l^SSl 
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SCMF simulations are hampered by the general limitations of the Edwards models 
- chains can cross each other. However, there are ways to introduce the dynamical 
effect of entanglements at least in cases where the equilibrium configuration space is 
not affected by topological constraints.!^^ 

IV. 3. 2 Field-Based Dynamics 

Dynamic Density Functional Theory (DDFT) 

Dynamic density functional theories (DDFT) are based on density functionals for 
polymer systems, such as, e.g., Eq. (j3.'25p . They supplement them by a model for 
their dynamical evolution at nonequilibrium. For diffusive dynamics, the dynamical 
equations in an imposed flow profile v have the general form.^ISSHIol] 

^ + V(vpi) = I dr'J^VrAiJ{r,r')Vr'J^ + r]^{r,t), (3.63) 

where Ajj(r, r') is a generalized mobility, and r]i(r, t) a Gaussian white noise with mean 
zero. If the amplitude of the latter is very small or zero, one has 'mean-field dynamics' 
and the system evolves towards a minimum of the free energy functional (although 
not necessarily the global minimum). If the noise is larger and satisfies the fiuctuation 
dissipation theorem, 

{ijiir,t)r],ir',t')) = -2k^T 5{t - t')V,Kij{v,v')V , (3.64) 

the density configurations {/9j(r)} in an equilibrium simulation (no flow, sufficiently 
long runs) will be distributed according to P[p\ oc exp(— J^[p]/A;gT). 

The kinetic Onsager coefficient Ajj(r, r') depends on the microscopic dynamics in 
the system. Since it characterizes the current of component i in response to an external 
force acting on component j, it is reasonable to assume that it is proportional to the 
local density Pi{r). An efficient choice which however disregards the connectivity of 
the chains is thus Ajj(r,r') = Mj Pi{v) 5{r — r') 6ij for compressible systems ,'2^ or 
A(r, r') = M paPb '5(r — r') for binary incompressible systems (local dynamics) 

To account for the chain connectivity, one must include information on the chain 
correlations. For example, for Rouse chains, AiAr, r') should be proportional to the 
pair correlators Kij defined in Eq. (I3.34p .^^^!^^ At first sight, using such complex 
Onsager coefficients in a simulation seems forbidding, but thanks to a clever trick due 
to Maurits and Fraaije,'^ it becomes feasible .'^'^ISlllI] 

The DDFT method has been extended, e.g., to account for viscoelasticP^S or hy- 
drodynamic eff ect s .IMllMIIIMSI ^ lattice version has also been proposed.'^^ 

Time Dependent Ginzburg-Landau (TDGL) Theories and Cell Dynamics 

Like dynamic density functional theories, time-dependent Ginzburg-Landau (TDGL) 
theories supplement a free energy functional of an 'order parameter' field <I>(r) by 
a model for the dynamical evolution of $. The TDGL theories of interest in multiphase 
polymer systems mostly operate with locally conserved order parameters and have the 
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same general structure than Eq. (j3.63p . For example, time-dependent Flory-Huggins- 
de Gennes theories are used to study the demixing dynamics in polymer blends, and 
time-dependent Ohta-Kawasaki theories to study the ordering kinetics in copolymer 
systems. Discrete lattice versions of TDGL theories are often referred to as 'cell dy- 
namics' models. 

We consider specifically the time-dependent Ohta-Kawasaki theory. Starting from 
the free energy functional (j3.42p for melts of diblock AB copolymers and choosing an 
Onsager coefficient that descri bes local diffusive dynamics, A(r,r') = M 6{r — r'), we 
obtain the dynamical equation^^SQI 



dt 



+ V(v$a) = MA 



Po 
N 



A 



r]{r,t) 

BA^a) - A{^A - ^a)] + r/(r, t), (3.65) 



where ^a is the total volume fraction of monomers A in the system. Hence the some- 
what awkward long-range 'Coulomb' term in Eq. (I3.42P becomes short range, and the 
dynamical equations only depend on local terms. Because of the appealingly simple 
structure of the final theory, Eq. (I3.65p . it is widely used for simulations of copolymer 
systems at equilibrium and under shear (see below) . Bahiana and Oono have formulated 
a discrete version on a lattic^^^ which is equally popular in cell dynamics simulations. 

TDGL simulations are much faster than molecular field simulations, but of course, 
the underlying model is less accurate. Honda and Kawakatsu have recently proposed 
a multiscale hybrid method that combines the two approaches, using dynamic density 
functional input to improve on the accuracy of the TDGL model.lSll Such hybrid 
approaches will presumably gain importance in the future. 



IV. 4 Applications 

After this overview over the main models used for simulations of multiphase polymer 
systems, I will now illustrate them by reviewing simulation work that has been done on 
homopolymer blends and copolymer melts. I focus on simulations of generic models. 
Atomistic studies or studies of bottom-up models are scarce and have already been 
discussed earlier (Sees. [rvTT] and irV.2.ip . 



IV. 4.1 Homopolymer Blends 
Bulk Properties 

I have already mentioned the pioneering simulations of binary blends by Sariban 
and Binder and by Cifra et a/.llSsHlHZl Following up on this work, a number of stud- 
ies, mainly by Binder and coworkers, have considered the critical behavior of binary 
blends. SS'^^^^^^^ As discussed in Sec. IHL2.5t fluctuations shift the critical demixing 
point and change the critical exponents from mean-field-like to Ising-like. However, the 
region where the critical behavior deviates from the mean-field prediction shrinks with 
increasing chain length, and effective mean-field behavior could be observed already for 
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relatively moderate chain lengthsP^ Coming from the other end, field-based Monte 
Carlo simulations of the Flory-Huggins-de Gennes functional ()3.13p have confirmed that 
noise shifts the coexistence curve and changes the critical behavior from mean-field to 
Ising.l^ 

Along with these studies of static blend properties, extensive work has also been 
dedicated to the dynamics of demixing. On the particle-based side, this was mainly 
investigated using kinetic Monte C ar lo .[2l2ll258ii259j j^g^g^er et al have compared results 
from kinetic Monte Carlo simulations and different versions of the dynamic density 
functional theory in a study of spinodal decomposition in symmetric blends.!^^ Other 
field-based simulations have mainly relied on time-dependent Ginzburg Landau models, 
which have the advantage that one can reach much later stages of demixing. They 
were used to study demixing processes at equilibrium^^^^^*^ ^^'^ and under shear PSHEMl 
Particularly interesting morphologies can be obtained if the dynamics in the two phases 
is distinctly different, i.e., one component becomes glassy during the demixing procesd^ 
or crystallizesPEOlEMl 

The particle-based studies mentioned so far have used coarse-grained models of 
blends that demix explicitly for energetic reasons. A number of authors have explored 
other factors that are believed to affect the chain miscibility with generic models,'^^ 
e.g., the effect of nonrandom mixing,'^^'^^ shape disparity,'^^ stiffness disparity,'^^ 
different architectures,'^^ and a different propensity towards crystallization.'^SlEzl 

Internal Interfaces 

In the miscibility gap, polymer blends have a highly inhomogeneous structure with 
droplets of one phase dispersed in the other phase. Their material properties are 
largely determined by the properties of the interfaces separating the two phases. While 
the distribution, size, and shape of the droplets depend on how the blend has been 
processed, the interfaces separating them often have time to reach local equilibrium 
and can be studied by means of equilibrium simulati ons . The first study of an interface 
in a binary blend was carried out by Reiter et alW^ in 1990. Since then, several 
authors have investigated interfaces in symmetricP^'^^HlEMEZZI or asymmetricP^SHHo] 

binary blends {e.g., blends with stiffness and/or monomer size disparities) by means of 
generic particle-based simulations. 

The local interfacial structure is of interest because it determines the mechanical 
stability of an interface. For example, the local interfacial width gives the volume 
in which chains belonging to different phases can entangle. On the other hand, we 
have already discussed in Sec. IIII.4I that interfaces exhibit capillary waves, which are 
significant on all length scales. This becomes apparent from the fact that the capillary- 
wave contribution to the total width, as given by Eq. ()3.53p . diverges both if the system 
size L becomes very large and if the microscopic coarse-graining length oq becomes very 
small. Therefore the length scales of the capillary waves and those of the local interfacial 
profiles cannot be separated clearly, and it is not clear, a priori, whether the concept of 
a 'local interfacial structure' is at all meaningful. This question has been investigated 
by Werner et ai^^® by simulations of the bond- fluctuation model (see Sec. IIV.2^ . They 
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demonstrated that it is indeed possible to describe homopolymer interfaces consistently 
in terms of a convolution of 'intrinsic' profiles with capillary wave undulations. They 
also studied the influence of confinement both on the capillary wave^^SUlgsai a,nd on the 
intrinsic interfacial width.'^SS 

Equilibrium simulations give information on the stability of interfaces under me- 
chanical stress, but in a rather indirect way. A few authors have probed directly the 
rheological properties of interfaces with nonequilibrium particle- and field-based sim- 
ulation methods, looking, e.g., at shear thinning and interfacial slip.'^^*"^^'' Detailed 
simulations of nonequilibrium interfaces are expensive, but with the development of 
new efficient simulation algorithms and modern fast computers, they become feasible. 



Surfaces and Films 

Another topic discussed intensely in the literature is the behavior of blends in the 
vicinity of surfaces, and that of blends confined to thin films. From a simulation point 
of view, these two situations are identical, because surfaces are typically studied in slab 
geometries {i.e., with periodic boundaries in two directions and free boundaries in the 
third). 

A large amount of work has been dedicated to the somewhat artificial situation of 
surfaces to which both types of monomers have exactly equal affinity. Even though 
these surfaces are perfectly neutral, one component will typically segregate to them: 
In incompatible blends, the minority component aggregates at the surface in order to 
reduce unfavorable contact^^SHHlSol contrast, the minority chains are removed from 
the surface in miscible blend^^^ . In blends of chains with different stiffness, the stiffer 
chains are pushed towards the surface, because they loose less entropy there.'^^'^^ 
For the same reason, linear chains aggregate at surfaces in blends of linear and star 
polymer s.'^^ 

Cavallo et al have systematically investigated the phase behavior of films confined 
between neutral walls as a function of the film thickness. If the film is made thinner, one 
observes a crossover from three-dimensional to two-dimensional Ising b ehavior j2251l296| 
Fluctuation effects in thin films are observed to be much stronger than in the bulk, 
consistent with our discussion in Sec. IIII.2.5t In two dimensions, the Ginzburg param- 
eter no longer scales with the chain length and stays finite for all chain lengths. A 
second transition occurs when the film becomes so thin that polymers are effectively 
two-dimensional, i.e., they can no longer pass each other. This reflects the theoret- 
ically expected fundamental difference between the demixing behavior of overlapping 
and non-overlapping two-dimensional polymers.'^ 

The more realistic situation of selective walls, to which one component adsorbs 
preferentially, has been addressed as well by different authors. In this case, the phase 
behavior is governed by wetting phenomena and capillary condensation.^^'' 

The studies discussed so far were based on particle simulations. A few authors 
have used field-based simulations to explore dynamical aspects of phase separation in 
thin polymer blends. Morita et al have studied the interplay of spinodal decomposition 
and interfacial roughening due to droplet formation with dynamic density functional 
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simulations Shang et al have used a time-dependent Ginzburg-Landau approach to 
study the spinodal phase separation of a thin film on a heterogeneous substrate.'^^ 

IV.4.2 Copolymer Systems 

Copolymers as Compatibilizers 

Copolymers were originally designed as natural surfactant molecules that increase 
the miscibility of incompatible homopolymers and enhance their interfacial properties. 
They are usually much more expensive than their respective homopolymers, but adding 
a small amount of copolymer can already improve the properties of the homopolymer 
blend significantly. 

A number of researchers have considered the effect of copolymers on the demixing 
transition for different copolymer architectures. ^^^ '^''^ Dadmun and Waldow^^^ have 
pointed out that copolymers not only shift the transition point towards higher values 
of the Flory-Huggins parameter but also change the critical exponents via a Fisher- 
renormalization mechanism. In the phase-segregated regime, copolymers aggregate to 
the interface, reduce the interfacial tension, and enlarge the interfacial width. The 
interfacial structure of homopolymer interfaces with adsorbed copolymers has been 
explored in detail by several authors .'^SIEMllsolHloi] 

Milner and Xi have noted in 1996 that the main compatibilizing effect of copolymers 
probably has a kinetic origin:'^''^ Copolymers reduce the rate of droplet coalescence dur- 
ing the processing of the blend via a Marangoni effect: If the copolymer concentration 
drops somewhere at the surface of a droplet, the surface tension increases locally. This 
induces flow in the direction of the weak point. Hence the copolymer film stabilizes 
itself kinetically, much like a soap film. In addition, the copolymer blocks stretching 
into the bulk form repulsive coronae. Experimental studies suggest that the resulting 
steric repulsion between droplets may be even more prohibitive for droplet coalescence 
than the Marangoni effect.l^^ Kim and Jo have studied the influence of copolymers 
on the dynamics of demixing in a series of works ,'3^2H3llI i^^^ they used kinetic Monte 
Carlo, hence they did not include hydrodynamic effects and could not study the effect 
of the Marangoni flow. 

If one increases the copolymer concentration beyond a certain threshold, the macro- 
phase separated phase at high eventually gives way to a microphase separated phase 
(see, e.g., the phase diagram shown in Fig. 13. 6p . A few authors have explored the full 
phase behavior of ternary systemsj ^^'^ ll ^^^ ^*^^^^^ In ternary systems containing A and C 
homopolymers and ABC triblock copolymers, a whole zoo of new tricontinuous gyroid 
phases can be observed. ^^"^^ 

Pure Bulk Copolymer Melts 

The propensity of copolymers to self-assemble into complex mesostructures makes 
them attractive for various micro- and nanotechnological applications, which is why 
pure copolymer melts have become interesting in their own right. Many simulation 
studies are now concerned with the properties of pure copolymer melts. 
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Figure 3.8: Phase diagrams for diblock copolymer melts as obtained from Monte Carlo 
simulations of a lattice model (chain length = 30) by Matsen et al. (leftp^ and 
from field-theoretical Complex-Langevin simulations (Ginzburg parameter C = 50) by 
Lennon et al (right). Symbols L,C,G correspond to diblock phases of Fig. 13. 3t in 
addition, the Monte Carlo phase diagram features a perforated lamellar (PL) phase. 
Left figure: Courtesy of Mark W. Matsen, adapted from Matsen, Griffiths, Wickham, 
and Vassihev, J. Chem. Phys. 124, 024904 (2006). Right figure: Reprinted with 
permission from Lennon, Katsov, , and Fredrickson, Phys. Rev. Lett. 101, 138302, 
2008. Copyright (2008) by the American Physical Society. 



Early studies of fluctuations and chain correlations near the order/disorder tran- 
sition (ODT), coming from the disordered sid^^HBEIHlETll j^ave reproduced the ODT- 
singularity in the structure factor and revealed the dumbbell structure of the chains 
mentioned earlier.HD Later, intense work has been devoted to studying ordered lamellar 
structures below the ODT in melts of symmetrical diblock copolymers, and analyz- 
ing them with respect to their structure, their dynamics, and their fluctuations .'^^QHSU] 
Simulation studies of asymmetric diblock copolymer melts have also reproduced most 
of the other mesophases - the cylindrical phase, the bcc sphere phase, and even the 

gyroid phase.204 211,235;M35] 

Locating the actual position of the ODT accurately is difficult, especially in lattice 
modelSjl^^ since the natural periodicity of the structures is in general incompatible with 
the box size. This results in complex finite size artefacts. Nevertheless, phase diagrams 
have been obtained in recent yearSj.'^^^^^^SSl Particle-based and field-based simulations 
have provided evidence that for symmetrical diblock copolymers (/ = 1/2), the transi- 
tion to the lamellar phase is shifted ,'2231l24QI compared to the mean-field phase diagram, 
and becomes first order ,'232H342j -^^ agreement with the theoretical expectatioii^^ (see 
Sec. lIII.2.51) . 

Two examples of phase diagrams obtained with different simulation methods are 
shown in Fig. l3.8t One was calculated by Matsen et a Pausing Monte Carlo simulations 
of copolymers with length = 30 in a simple lattice model (left), and one by Lennon et 
al¥^ using field-theoretical calculations (right) at Ginzburg parameter C = 50. Both 
phase diagrams significantly improve on the SCF phase diagram of Fig. 13.51 (right) 
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in the weak segregation regime, and reproduce the main quaUtative features of the 
experimental phase diagram (Fig. 13.51 left): The transitions are first order everywhere. 
The ODT is shifted to higher x^- Direct phase transition between the disordered phase 
and the complex mesophase (PL or G, respectively) or the lamellar phase are possible 
for a range of copolymer block fractions /. The phase diagrams even reproduce a small 
hump of the ODT transition line at the boundary to the complex mesophase (PL or 
G), which is also observed experimentally. The only 'problem' with the Monte Carlo 
phase diagram is that it features a perforated lamellar (PL) phase instead of a gyroid 
(G) phase. This may be a finite size artefact, as suggested by the authors, or a property 
of the lattice model under consideration. (Gyroid phases have been found in lattice 
modelsjl^^ but it should be noted that the free energies of the PL and the G phase are 
very close according to SCF calculations.) Nevertheless, the simulations demonstrate 
convincingly that the discrepancies between the experimental phase diagram and the 
SCF phase diagram shown in Fig. 13.51 can to a large extent be attributed to the effect 
of fluctuations. 

In recent years, researchers have also begun to simulate melts of more complex 
copolymers, e.g., starblock copolymers,'^^^^^^^ rod-coil copolymers ,'^^'^^diblock copoly- 
mers with one crystallizing component [2lH26l|350H352j ^j-ii^iQ^k cop olymer s I2521{356] xg^n- 
dom cop oly mer s .I^SIIISSH] interesting study on diblock copolymers with one am- 
phiphilic block has recently been presented by Khokholov and Khalatur.l^^ Since the 
amphiphilic block favors interfaces, the morphology of the mesophases changes com- 
pletely and is characterized by thin channels and slits. Lay et al. and Palmer et a/ have 
studied the computationally challenging problem of microphase separation in randomly 
crosslinked binary blends, '^^'^ '^^^ and also the inverse problem, to which extent ordered 
copolymer structures can be stabilized by cross-linking.'^^ 

A large amount of simulation work on copolymer melts has been done with time- 
dependent Ginzburg-Landau approaches j^MlEzSI xhese studies have mostly addressed 

dynamical questions, i.e., the kinetics of ordering, disordering processes in pure 
and in mixtures containing copolymers i^ZSlSzHl Already in pure diblock copolymer melts, 
ordering/disordering processes were found to proceed via intricate pathways that in- 
volve nontrivial intermediate states {e.g., the perforated lamellar state which is only 
metastable at equilibrium). In mixtures, the interplay of macrophase and microphase 
separation leads to a wealth of new transient morphologie^^^^^^^^ (see also Ref. i379|) . 

Confined Copolymer Melts 

In recent years, there has been growing interest in confined copolymer systems. 
The first studies have explored the ordering of copolymers melts in thin films between 
neutral wallP^^^^^^ or general (selective) walls ,'^^^'^294] both with particle-based models 
and dynamic density field theories. Triblock cop olymer ^^^HSoolaj^d copolymer blend^^ 
have also been considered. The dynamics of copolymer ordering in confinement was 
studied with time-dependent Ginzburg-Landau methods )2921I4Q1J|402] ^ series of papers 
have dealt with the technologically relevant question, whether and how surface patterns 
can be transferred into copolymer films, .^^S9|| 403H408 | 
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The current focus shifts to nanocyhndricany or sphericahy confined blends P^S^SHUS 
A variety of new structures can already be observed in systems of diblock copolymers 
confined to nanocylinders, e.g., mesh structures, single and double helices.'^lSl For tri- 
block copolymers, the spectrum is even more diverse.!^^ The confined structures depend 
on the bulk structure and on the shape of the confining channels. Li and coworkers 
have shown that possible morphologies can be screened efficiently with the method of 
simulated annealing.SHlEl] 

Copolymers under Shear 

Finally in this section, I briefly mention the large amount of work that has addressed 
the effect of shear on the microstructure of copolymer systems, with both field-based 
and particle-based methods I^^SHMU Special attention was given to the phenomenon that 
lamellae reorient under shear .ll2ll|442H445| Experimentally, it is observed that lamellae 
orient parallel to the shear flow at low shear rates, and perpendicular to the shear flow 
at high shear rates. Simulations give conflicting results. The parallel state consistently 
becomes unstable at high shear rates, but its relative stability at low shear rates (or 
results on the indicators of r elative stability such as the entropy production) seems to 
depend on the type of model.HSllMSj 

V Future Challenges 

The equilibrium theory of fluid polymer mixtures is fairly advanced. Thanks to the uni- 
versal properties of polymers, it requires relatively little input information (x parameter 
etc.) to be predictive at a quantitative level. However, if one goes beyond equilibrium 
and beyond the fluid state, the situation is much less satisfying. For example, very 
little work has been done on crosslinked polymer blends ,'^^2'^^^ even though they are 
common in applications. Physical crosslinks can be established if domains crystallize 
or become glassy. As we have seen above, research on blends with one crystallizing or 
glassy component is also rather scarce. 

Another important issue is the influence of the blend processing on the proper- 
ties of the resulting materials, i.e., the structure of phase separating blends under 
shear. The mechanical properties of immiscible blends depend crucially on their mi- 
croscopic morphologies, i.e., the sizes and shapes of droplets, which in turn depend 
on the manufacturing process. Theoretical nonequilibrium state diagrams that relate 
the processing conditions (shear rates, geometry, copolymer concentration etc.) with 
final morphologies are still missing. Since the relevant length scales are relatively 
large and hydrodynamics are important, the simulation method of choice should be 
a cell dynamics method that incorporates hydrodynamics, e.g., a Lattice-Boltzmann 
method.'^^ Methods that combine Ginzburg-Landau functionals for immiscible fluids 
with Lattice-Boltzmann models for Newtonian fluids have been developecP^ and used 
to study demixing processes at reslpSl and under shear .S^SHiSU Qij-^ud and cowork- 
ers have proposed a Lattice-Boltzmann method for viscoelastic fluids, which is more 
suitable to describe polymers, 12221 and carried out first simulations of viscoelastic liquid 
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mixtures Nevertheless, the whole field is still in its infancy. 

As the field of polymer simulations reaches maturity, the bottom-up modeling ap- 
proach (Sec. lIV.ip will gain importance. So far, the vast majority of theoretical and 
simulation studies of (co)polymer blends was based on generic model systems. One or 
two decades from now, realistic simulations of specific polymer blends will probably 
be equally, if not more common. One of the major challenges in this context is to 
develop hybrid multiscale methods that combine different levels of coarse-graining, i.e., 
use a relatively coarse basic model and fine-grain selectively in interesting regions of 
the materials {e.g., interfaces) 
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